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Abstract 

Theoretical and experimental investigations of upper atmosphere dynamics 
have been pursued under this grant for almost two decades. This final re- 
port presents a brief overview of the significant contributions made to the 
understanding of the dynamics of the earth's upper atmosphere, including 
the addition of winds and diffusion to the semi -empirical Global Reference 
Atmospheric Model developed for the design phase of the Space Shuttle, 
reviews of turbulence in the lower thermosphere, the dynamics of the equa- 
torial mesopause, and stratospheric warming effects on mesopause level 
dynamics, and the relevance of these studies to the proposed Middle Atmo- 
sphere Program (1982-85). The report ends with a chronological bibliography, 
including abstracts, of all papers published under this grant. 
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CHAPTER I 


Overview 

In the early sixties, the foremost technique for the measurement of 
wind profiles in the lower thermosphere was the rocket released chemical 
vapor trail. The early work performed under this grant consisted of the 
development of techniques for the accurate measurement of trail position 
with time (based on triangulation against the stellar background from two 
or more ground based camera sites), and lead to the determination of wind 
and wind shear profiles, gravity wave and turbulence spectra, and represent 
ative turbopause heights. In the mid and late sixties, better models ex- 
plaining the turbulence observed evolved, supported by concurrent interpre- 
tation of wind shears in the 80 to 100 kilometer height range, measured by 
the multi -station meteor wind radar technique. The cascade of energy from 
the diurnal tide to gravity waves to turbulence was proposed at this time, 
and subsequently verified by others. 

Work on turbulence in the lower thermosphere continued through the 
early seventies, culminating in the review presented here as Chapter III. 

Also during the early seventies, the space shuttle reentry model atmo- 
shpere, and its extension to a wind and diffusion model, were developed. 
This Global Reference Atmospheric Model (GRAM) is available on magnetic 
tape, and has been widely used both within and outside NASA. The model 
describes pressure, density, temperature and winds as monthly means and 
variances (details in Chapter II). 

Funds from this grant have been used in support of the analysis and 
interpretation of winds from the Georgia Tech Radio Meteor Wind Facility. 
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This is the only continuously operating meteor wind radar In the world which 
produces height/time profiles of winds In the 80 to 100 km height range (a 
sample of the prevailing winds, zonal and meridional, for the period August 
'74 - August '75 Is presented as Figure 6 In Chapter VII, paper b). 

Mesospheric and lower thermospheric dynamics have been Investigated 
(Chapter IV), and more recently. In cooperation with the French National 
Center for Telecommunications Studies (CNET), the dynamics of the equatorial 
mesopause (Chapter V). 

Of particular interest to meteorologists and aeronomers In recent years 
has been the delineation of the coupling between various height ranges 
within the atmosphere. One obvious source of this coupling Is to be found 
In random internal atmospheric gravity waves. Their contribution to vertical 
mixing In the stratosphere and mesosphere has been Investigated, and Is de- 
tailed In Chapter VI. 

More recently, the effects of winter polar stratospheric warmings on 
the circulation at midlatitude mesopause levels has been studied, and pre- 
liminary results of this Investigation are presented In Chapter VII. 

This report would be incomplete without a reference to MAP - the Middle 
Atmosphere Program, scheduled for 1982-85, which Is discussed In Chapter VIII. 

Finally, and probably most Importantly, a bibliography of all publica- 
tions produced under this grant, together with abstracts. Is appended. 


CHAPTER II 

The Global Reference Atmospheric Model 

Although the Global Reference Atmospheric Model (GRAM) was developed 
as a thermodynamic variable model irt support of the design phase of the 
Space Shuttle by a grant from NASA Huntsville, part of its subsequent evo- 
lution into a wind and diffusion modsl was supported under this grant. 

GRAM is a computer program, available on magnetic tape, which combines 
the 4-D model of Spiegler and Fowler (1972) (0-25 km), a development of the 
Graves (1971) model (for the region between 30 and 90 km), and the Jacchia 
(1970) model (115-700 km). Smooth transitions between these models are 
accomplished with a fairing technique (90-115 km) and an interpolation 
scheme (25-30 km). 

In addition to monthly mean values of pressure, density, temperature 
and winds, two types of perturbation are evaluated; quasi -biennial (QBO) 
and random. For further details of this model, see "A Global Reference 
Atmospheric Model for Surface to Orbital Altitudes" by C. G. Justus, R. G. 
Roper, Arthur Woodrum and 0. E. Smith, Journal of Applied Meteorology, 1^, 

3-9, 1976. 
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Turbulence in the Lower Thermosphere 


(Chapter 7 of "The Upper Atmosphere and Magnetosphere," published by 
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Turbulence in the 
Lower Thermosphere 


ROBERT C. RORER 

Georgia Institute of Technology 


7.1 PROLOGUE 

A static, or motionless, atmosphere can be adequately 
described on a macroscopic scale by three parameters; 
pressure, temperature, and composition. If the atmo- 
sphere were macroscopically motionless for a sufficiently 
long time (its molecules in purely random motion), one 
would (ind a decrease in mean molecular weight with 
altibide, the lighter species having separated from the 
heavier species by the process of molecular diffusion. 
Such a process is not observed in the earth’s atmosphere 
until one reaches an altitude of from 1(N) to 1 10 km; ap to 
this altitude the mean molecular weight remains approx- 
imately constant. This comes about because the earth’s 
atmosphere is not motionless on a macroscopic scale but 
is continually mixef* by atmospheric turbulence and 
large-„cale circulation. The beginning of diffusive separa- 
tion at 100-1 10 km does not mean that there are no winds 
to cause turbulence and mixing above these altitudes but 
rather that at such low densities molecular diffusion rates 
are much greater than the wind-induced mixing rates. 


this chapter concentrates or the region known as the 
lower thermosphere, between the mesopause at approxi- 
mately 80 km, where the atmosphere is always mixed, and 
the altitude of 130 km, where diffusive separation is 
always observed. 

Cur knowledge of mixing processes in the upper at- 
mosphere has been accumulated mostly in the years since 
World War 1*. While the lower thermosphere has been 
probed for decades using radio techniques, detailed 
knowledge of itr structure has only come with the use of 
rocket soundings. 

Before considering the techniques of measurement of 
mixing processes in ^e upper atmosphere and the inter- 
pretation of those measurements, a deffnition of "mixing 
process" is appropriate. On a global scale, the atmosphere 
can be said to be mixed by the transport of constibients 
from one location to another by large-scale wind systems. 
This mixing is important. Here, however, we will concern 
ourselves with more localized mixing, usually brought 
about in the free atmosphere by shears in the wind system 
and called atmospheric turbulence. A turbulent atmo- 
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iphorv is miied at rulL's often hundreds ur more Umes 
faster than its muleeules van diffuse by means of their 
thermal metiuns, and thus hirbulenec ean be very eflee* 
hve in maintaining homiigeneity in an atmosphere lom- 
pused of several molecular speeiesj h< nee the term 
“homosphere” is often applied to deseribe the homoge* 
neous atmosphere, with the heterosphere above. The 
level at Whieh turbulent mising eeases to be effective 
is sometimes referred to as the homopause, but here we 
will use the term turbopause 
In addition l.'> mixing tlie atmosphere, locally intensive 
turbulence also eouses local heating. Turbulence is dis- 
sipative, extracting energy from Uie total flow and brans* 
ferring it by a cascading process to scales so small that the 
random motion of the moleiules, which determines the 
temperature, is increased. Thus turbulence in the lower 
tlier,T)osphere is important because its intensity affects 
both the relative concentrations of constituents of the 
thermosphere and the heat budget in the lower thermo- 
sphere. The source of the turbulent energy resides ulb* 
mately in the lower atmosphere, and this energy is trans- 
ferred to the thermosphere by the upward propagation of 
internal atmospheric gravity waves and ti^l winds (see 
Chapters 8 and 0). 


7,2 INTBODUCTION 

In the parlance of ionospheric physicists, the lower ther- 
mosphere is known as the E region. The E region has 
been probed from the ground almost since the inception 
of radio. In particular, irregularities in the structure of the 
ionization that produces reflections of radio waves at this 
altitude have been observed for decades. Because these 
irregularities could be observed only through reflection 
fronri the continually changing ambient ionization imbed- 
ded in the neutral atmosphere, their interpretation in 
terms of the turbulent structuri' of the neutral atmosphere 
was highly speculative. However, in the early 19^'s, a 
technique was devised that enabled the characteristics of 
the neutral atmosphere to be determined by direct obser- 
vation of relatively well understood radio reflectors — the 
ionized trails formed by meteorites entering the earth's 
atmosphere. 

Measurements of the radio-frequency Doppler shifts 
produced as meteor trails were blown along by the neu- 
tral winds between 80 and 100 km (the altitude range over 
which must meteorites bum up in the earth's atmosphere) 
were, and still are, used to provide knowledge of the 
neutral wind and its variation in botli height and time. At 
the same time, considerable progress was being made by 
scientists working in the analysis and interpretation of 
atmospheric turbulence in the troposphere. The interpre- 
tation was assisted in no small measure by the progress 
made in turbulence theory, in particular by Batchelor's 
(1053) Theory of Homogeneous Turbulence. After this 
theory was applied successfully to the explanation of the 
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tropospheric scattering of veo’-high-frequency radii, 
waves, Booker and Cohen (1056) attempted to explain the 
fading observed on lung-duration meteor echoes in terms 
of turbulence in the neutral atmosphere at E-region al- 
titudes. From their data, tliey deduced tiiat energy was 
extracted from tlie large-scale wind motions at meteor 
altitudes and dissipated at a rate s - 25 \k'/kg. While the 
underlying theory was sound, their paper was attacked on 
the basis of their interpretation of tlie echo-fading pro- 
cess. 

In the late 105O's, the chemical-release racket tech- 
nique was perfected and used to introduce a visible 
tracer, initially sodium, into the atmosphere over the 
altitude range 80 to 200 km. Such a release, made at 
twilight white the trail was illuminated by sunlight and 
the ground was in darkness, could be photographed from 
several camera sites on the ground, and a time series of 
exposures recorded simultaneously could be used for 
triaivgulation of tlie release, thus detennining its motion 
with time. Winds could be determined fur as long as tin- 
trail remained visible, sometimes for as long as IS min. 
Since the early 1960's, trimethyl aluminum (tma), which 
reacts with atomic oxygen in the ambient atmosphere to 
produce products that not only scatter sunlight but also 
glow in the dark (making nighttime as well as twilight 
measurements possible), has been used in addition to 
sodium. For daylight releases in particular, lithium has 
been used. 

One outstanding feature of these trails is the fact that 
below an altitude of some 105 km tlie release is "obvi- 
ously turbulent," whereas above that altitude the trail 
expands smoothly, as it should under the action of 
molecular diffusion alone. Considerable controversy still 
exists over the interpretation of data from rocket-released 
vapor trails as evidence for turbulence of the ambient 
atmosphere. Some interpretations of sodium vapor trails, 
for example, have led to anomalous results that are 
thought to be related to the energetics of the thennite 
bum used to produce the sodium vapor. However, in- 
terpretation of the breakup of vapor trails below what has 
become known as the turbopause is not the only evidence 
for the existence of turbulence at tlicse altitudes. Theoret- 
ical studies of the atmospheric heat budget in the high 
atmosphere by Johnson and Gottlieb (1070), for example, 
require vertical diffusivlties below tlie turbopause that 
are much larger than molecular. Calculations of the dif- 
fusivities responsible for the measured constituents 
above the turbopause also require similar values of dif- 
fusivity below the turbopause; composition mea- 
surements made using rockelbomc mass spectrometers 
show the level of diffusive separation to be considerably 
higher than would be the case in a nunturbulent atmo- 
sphere. Further, a completely different class of mea- 
surements, based on the shearing of radio meteor trails by 
the winds in the 80- to 100-km region ratlier than the 
fading of individual meteor echoes, yields values consis- 
tent with those deduced from vapor-trail observations. 
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7.3 THE INTERPBETATION OF 
MEASUREMENTS 

li it pnssiblf fur rm'ki>t-«vukf i-fR'i'tii and vheinlcal energy 
releaned by enniuniiimnik (from the thennile canister) to 
innuente tlie snbseiiuent dispersion? Certainly. Such 
efTects have been documented. For energetic releases 
such as rwliet eshaiists or large quantities of esplosivc, 
tlie "release phase" at thermospheric altitudes lasts only 
fur some 10 sec or less. Figures 7.1 and 7.2 show an 
escellent example of a release phase anomaly in a TMa 
trail, as photographed by the Smithsonian Institution’s 
Baker'Nunn camera at Woomera, Australia (31* S). Figure 
7.1 shows 0 portion of the trail 2 sec after release in the 
100. to 1 lO-km height Vangei vortex shedding on the 



FIGURE 7.1 Portion of a brimelhyt aluminum trail showing a 
release phase instability with vortex shedding on the right-hand 
edge 2 see after release at 103 km. 



FIGURE 7.2 Same Os for Figure 7.1. but 6 see after release. 
Note that the release phase instability has been damped out 


right-hand side of the trail is clearly visible. In Figure 7.2, 
a frame taken 4 sec later, the efTects of this motion have 
disappeared. The subsequent breakup of the trail, with 
the production of the characteristic "obviously turbulent” 
appearance, did not occur until some SO sec later, os 
illustrated in the isodensitrace montages of Figure 7.3. 
These unique Baker-Nunn photographs are an example of 
the rewards of international cooperation. The Smithso- 
nian Institution for many years operated a worldwide 
network of Baker-Nunn cameras for the photographic de- 
termination of satellite positions. Relationships between 
the satellite station staiT and the Australian and British 
rocket experimenters at Woomera were such that the 
Baker-Nunn camera would be used to photograph the 
rocket releases when such use did not interfere with the 
primary mission of the observatory. 


original 
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FICUnE 7.3 Time history of a poitiun of the trail of Figure 7.2 showing the initial laminar behavior and subsequent 
breakup indueed by ambient atmospherie turbulenee. Quantitative measurements of dispersion ean be made from the 
isodensitraee eontours. The kink is at an altitude of 102.7 km. Note that the trail is laminar at +30 see and that the trail 
diameter at this time is greater than the scale of the eddies at +S0 sec. 


The abrupt cutofT in ambient turbulence (the tur- 
bopau.sc), which almost always underlies a region of high 
wind shear, is illustrated in Figure 7.4. That the bre^* 
down from "obviously laminar" to "obviously turbulent" 
represents a dramatic change, easily recognized by visual 
inspection of the film, is shown in Figure 7.S, in which 
the srpiare of the eflective radius of the trail is plotted 
against time after release at 105 km. The growth in the 
first few seconds after release is an order of magnitude 
faster tlian the subsequent growth up to the time of trail 
breakup, and this is regarded as representing the release 
phase, in which the energetics are deftnitely nonambient. 
Ihe growth between 8 and 32 sec after release could be 
molecular diffusion of a cloud with an inib'al radius of 130 
m, i.e., the release produced a cloud with, effectively, this 


radius at zero time. However, one can only say tliat the 
growth during this phase could be molecular — the 15 
percent error In the determination of the effective radius 
from each film frame, and the fact that the trail cross 
section is only approximately Gaussian, precludes a mea- 
surement of the diffusion coefticient to better than a factor 
of 3, The spreading due to turbulence commences 33 sec 
after release and then proceeds to follow tlie dispersion 
relation predicted by the theory of homogeneous turbu- 
lence; 

r,* a f », 

until, at 54 sec, the trait becomes too irregular for an 
estimate of radius based on a Gaussian distribution to be 
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Polnii A and Bo'oai 106 <nd 108 km altitude respeciMy, 0 1 3 3 km. 

FIGURE T.4 Montagi* nftliubrhaviuruflhu trail above and below the lurbopausv. Pointr A «nd B are at 106 and 108 km. reipecUvely. 
Note the tramition from laminar to turbulent at A between +23 and +61 see. Tliit has not oeeurred at B. 


tneanitiRfuI. However, Figure 7.5 does provide two pa- 
rameters: die effeetive radius at transition and its time of 
oeeurrenee after rcleafie, which should characterize the 
small-scale end of the turbulence spectrum. The concept 
is that until die trail spreads to a size equal to that of the 
smalle.st scale eddies, it is distorted but not spread by the 
eddies. Thus die scale size at transition is interpreted as 
indicating die size of the snialleBl eddies. 

Tlie most important parameter of any turbulence spec- 
tnim is e, die rate of dissipation of turbulent energy. In 
order to calculate c from any set of space-time correla- 
tions, it is necessar>‘tousea model.The simplest is that of 
A. Kolmogoroff, as elaborated on by Batchelor. Kcl- 
mugorofl' put forward die hypodiesis that, in any turbulent 
(low Reid in which energy is entracled from the large- 
scale or mean motion and cascaded to smaller scales 
before eventual dissipation at scales where molecular 
viscosity becomes important, there exists a range of scales 
siilllciendy removed from the large-scale anisotropic ed- 
dies, and yet not appreciably damped by molecular vis- 
cosity, that is both homogeneous and isotropic. The as- 
sumptions of homogeneity and isotropy inherent In this 
model are open to question in the case of lower thermo- 
spheric turbulence, but a considerable amount of work by 
several experimenters has shown diese assumptions to be 
reasonable, at least for length scales less than 1 km. 

Batchelor defines the basic tengdi parameter of the 
viscous region (that is, the size of the small-scale eddies 


that ore responsible for the conversion of the eddy kinetic 
energy to heat) as 

q - (>•»/€)«'•, 

where e is the kinematic viscosity (a measure of molecular 
dilTusivity), and q is the scale size la meters.'radian (dius 
the characteristic wavelength of these turbulent eddies is 



FIGURE 7.S Variation of (he ftqitare of Ute radius of (he trail of 
Figure 7,1 vkith time after release at 105 km. 
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FIGURE 7.6 Ptofilet of (he rat« of dittipaUan of lurbuicnt 
vrieigy deduced fium wind ahcao measured for 2S traits over 
Walliips Island (38* N) by Ruper (1966b) and from the dispersion 
of* similar number of trails over Eglin Air Force Base (2^ N) by 
Justus (1067), The upper scale gives a measure of the atmo- 
spheric heaUng resulting from tuniulent dissipation. 

2tt)). Batchelor also derines a t haracleristic time constant 
»• corresponding to diis unit length scale such that 

1) » (pr*)">. 

!• is the lifetime of the smallest eddies, essentially the 
time required for their dissipation by molecular viscosity. 
Combining these two eqiiations to eliminate i) yields 

»• « 

In tenns of the length scale r/, 

* a If*. 

In terms of the time scale I*. 

« or »•-*. 

Wind shears have been used to calculate the variation 
with height of the rate of dissipation of turbulent energy 
(Roper, 1966b). The kinematic viscosity v may be deter- 
mined from the viscusity and density published in the 


U.S. Standard Atmosphere Supplements for 1066. The 
frinctinn 

p (m* see*') ■ cap 10,17 (a • 80.0)1, 

with z in kilometers, fits the data at die tabulated al- 
titudes. The * profile determined in this way is shown in 
Figure 7.6, together with the profile inferred from mea- 
surements made by Justus (1967), calculated from the 
velocity fluctuations observed on 18 TSW trails. These 
values of < and e have been coupled to produce what can 
be regarded as average height-dependent characteristic 
length scales t)j| (• 2 mj, in order to express the length 
scale in the more usual wavelength notation, meters/ 
cycle). r)« and f* are platted, together with s'alues calcu- 
lated from Justus’s f profile, in Figures 7.7 and 7.8, 
respectively, Also plotted are the values determined from 
two Skylark-released TMA trails photographed by die 
Baker-Nunn camera at W’oomera at dawn and dusk on 
May 31, 1968. 

The similarity in form, and the order of agreement, 
between the predicted "average” trail time constants t«* 
and tf (Roper's and Justus’s values, respectively) and the 
measured values for each of the two releases is surpris- 
ingly good when one considers the approximations in- 
volved in (a) the model atmosphere viscosity, which is 
based on en average atmospheric model; (b) the variation, 
both diurnal and seasonal, in the turbulent dissipation 
rate, which has been averaged out in the construction of 
the t«* and t/ profiles; and (c) the fact that, while the 



FIGURE 7.7 The variation with alhhide of the KfllmogorolT 
miitoscales r/j and i)« inferred fium the turbulent dissipation 
proRles of Figure 7,6, Also shown are the microscales infeired 
from the May 1968 trails, r,* is the Gaussian radius of the 
evening trail at the time of onset of trail breakup, t*. 
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FiGl’RCTB Thi' vanatjnn wllhallitudeoflhe lime constant of 
the KciImpgatofT microscale for the f profiles of Figure 7.6 and 
the May 10^ trails. 

parameter t* appears to be a measure of the time tiikcn for 
trail breakup to oecut) It is not clear why it should be- 

The definition uf rj is a matliematical expedient charac- 
terizing tlie scale size for turbulence at which viscous 
dissipation becomes important, There is no obvious rea- 
son v,'hv uitlicr 7) or I * should be physically measurable 
featUY >' of the motion. Nevertheless, the close similarity 
in the sh.apes of the f* curves strongly indicates that the 
time delay in the unset of turbulence should be related to 
the time constant Of tlie Kolmogorolf microscole. Fur- 
thermore, tliese observations suggest a reason why tire 
lurbopause, dcfiited as the boundary between tlie regions 
that break up and those that remain laminar, should 
manife.st itself so abruptly. Above lOS km. tlie time con- 
stant for the Onset of turbulence increases so rapidly 
with altitude that the trail is not, in general, observed fora 
suflreient length of time for visible breakup to occur. 

Attempts have been made to explain the existence of 
the turbopause in tenns of a critical value of some param- 
eter, generally the Reynolds number R, or the Richardson 
number B|. These attempts have met with marginal suc- 
cess at best, partly because of the difficulty in defining 
the characteristic lengths tliat occur in these parameters 
and partly because it is not evident a priori what critical 
value the parameter should have at the turbopause. With- 
out a detailed knowledge of the temperature gradient at 
the scale characteristic of the vertical mixing process (a 
few hundred meters or less), the Richardson number just 
cannot be specified. Johnson (1975) has considered the 
relative importance of buoyancy and dissipation in some 
detail, referring particularly to the work of J. D. Woods, 
who determined that there was hysteresis in the 
criterion— laminar flows become turbulent when R|S 


0.25, while turbulent (lows become laminar when Ri si. 
While tlie mean (undisturbed) atmospheric temperature 
profile is pertinent to the breakdown from laminar to 
turbulent flow, once turbulence is established its cessa- 
tion will depend on tlie temperature grarient as modified 
by the torbulencc. As yet, there is no technique available 
fur measuring such temperature gradients in the lower 
thermosphere- 

Another enor that has often been made in attempting to 
explain turbulence in tie- lower tliennosphere is the as- 
sumption tliat the turhopause corresponds to an altitude 
at which turbulence ceases abruptly. The results pre- 
sented here, on the contrary, show that the turbopause Is 
the altitude at which the time constant of the Kohnugorun' 
microscale of the turbulence increases very rapidly with 
altitude. This viewpoint resolves the paradox that regions 
above the turbopause, which were tliought of as nontur- 
bulent, have difiusion coeflicients based on the measured 
laminar trail growtli that are greater tliaii molecular. We 
now sec that turbulence docs exist above the turbopause 
but that its efficacy in transport relative to molecular 
diffusion decreases with altitude. At an altitude of 130 
km, tlie contribution of turbulence to diffusivity is insig- 
nificant, even though its absolute value may be as large os 
it is at the turbopause, 

The rate of dissipation of turbulent energy may also be 
calculated from r), the length scale at which eddy diffu- 
sion becomes effective. However, since this lengtli de- 
pends crib'cally on the chape of the cloud (the assumed 
Gaussian variation across die cloud is rarely realized in 
practice below the turbopause) and because i oc r)-*, Rees 
et al. (1972) chose to use the relatively precisely deter- 
mined (* values to calculate the variation of < witli height. 
This is shown for a pair of dawn and dusk releases above 
Woomera t30" S) on May 31, 1968, in Figure 7.9. Note that 



FIGURE 7.6 Estimalc-s of tin rate of dissipation of turbulent 
energy • from the time of onsei of turbulence fur the release of 
May 1668. 
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du' liirliopaiisi- iii hit!hi‘r in the ewning (hun in liif 
niurniiig and llial Hit higln-r turbti|)auiif is assuiiatfd with 
a gri'alir iivt'rall tiirtiuU'nt inlt‘nsli>'. Utis siibstanlialfs 
tbr lugefslinn uf a variation in midlatitudo lurbupausi* 
hvigbi madf b>' Elford and itfipvr (1(167), uhiib was 
based on seasonal variations in turbulent intensit>’ at 63 
km as determined from the wind shear measured simul- 
taneously on individual radio meteor trails. 

The results from l\so further TslA releases made above 
Woomera at dawn on Ot tober 16 and dusk on Oetober 17, 
1666, are presented in Figure 7.10, For this pair of trails, 
the turbulent intensity is higher in the morning than in 
die evening— opposite to the May 1968 releases. This 
diurnal variation with season is die same as diat measured 
for the large-seale turbulent vripeit)' component from 
radio meteor winds at Adelaide (35* S). These Ottober 
releases are of particular interest in that they show alter- 
nating laminar and turbulent regions similar to those 
previously reported by Blamont.and Barat (1967). The 
various layers observed in these releases do not seem to 
be quite so simply related to the wind proflie as those of 
Blamont and Barat. However, die regions of prolonged 
laminar behavior all seem to be located at altitudes where 
the wind shear is high. 

In an attempt to explain why turbulence in the lower 
thermosphere should be stratiBed at times, and why in 
fait an ostensibly highly stable region of the earth's at- 
mosphere should be turbulent at all, Lloyd et al. (1672) 
proposed a model in which random internal gravity waves 
produce turbulence accompanied by a considerable mod- 
ification of the temperature profile. The creation of tem- 
perature inversions by turbulence is commonplace in the 
troposphere. It is proposed that a similar effect occurs in 



FIGUItE 7.10 Same as for Figure 7.0 but for the releases of 
October 1969. Note the presence of turbulent "sheets” alternat- 
ing with laminar layers. Diumally averaged dissipation rates for 
the month of October 1961 (Elford and Roper, 1967} and October 
1969 (McAvoney, 1970) as measured by the radio-meteor tech- 
nique at Adelaide (35* S) are shown for comparison. 


the mtm.* stable lower thermosphere, with part of the 
vertical ctimponent of the gr,-ivity-wave velocity ns the 
source of the liirhulcnt energy, Ifrandom, internal gravity- 
waves propagate at a significant angle to the vertical, then 
dieir vertical velocity component will be able to contri- 
bute to the destabilization of the stably stratified lower 
thermosphere (Hodges, 1667). Similar destabilization 
mechanisms have been discussed for the oceans by Phil- 
lips ( 1671 ) and for the lower stratosphere by Roach (1670) 
Radio meteor studies have already established that grav- 
ity waves are the source of die turbulent energy in die 
lower diermosphere. This has been confirmed quite in- 
dependently in radio-meteor studies by Sptzzichino 
(1672) and by Blamont and Barat (1667) through observa- 
tions of chemical releases. 

The lower thermosphere is stabilized against vertical 
motion by the mean temperature gradient: above the 
mesopause near 80 km the mean temperature increases 
with height, By equating the work done in the vertical 
displacement of a parcel of gas to the fraction of the 
vertical component of the gravity-wave spectrum that is 
responsible for the measured turbulence spectrum, Lloyd 
el ol. deduced the modification of the mean temperature 
profile that would be produced by the turbulence mea- 
sured on the October 16 trail in the height range 104 to 
110 km An Isodensitrace montage of this portion of die 
trail is shown in Figure 7.11. Tlie choice of the turbulent 
layers within which the analysis can be applied is some- 
what subjective, being based on regions where growth is 
"obviously" different from that above and below. 

The solid lines of Figure 7.12 show the results of the 
application of the Lloyd et al. model to each of the layers 
delineated in Figure 7. 1 1. The dashed lines are necessary 
for profile continuity and must represent laminar sheets. 
Because of the discontinuous nature of the determination 
of the modified profile and die forced fitting of the mid- 
point temperatures, the magnitudes of the positive and 
negative gradients are open to question. However, it is 
interesting to note diat the existence of similar gradients 
in the lower stratosphere, a region of similar mean- 
temperature gradient, is well documented from aircraft 
observations, as can be seen from the project hicaT de- 
termination shown in Figure 7,13 from Mitchell and Pro- 
phet (1969). Unfortunately, the flight of an aircraft at 
constant altitude cannot reveal a height profile of turbu- 
lence, but at least in the encounter with cat (clear air 
turbulence) at 21 km, the temperature profile has been 
modified in a manner consistent with the present model. 

Several deficiencies exist in the model, since the finite 
time constants of the processes involved (the period of the 
destabilizing gravity wave, for example) have not been 
considered. The basic energy-budget equation can be 
made more general by the inclusion of terms describing 
energy sinks (e.g.. heat conduction). One promising 
model being developed uses the reversible heating as- 
sociated with propagating gravity waves (Hines, 1965) as 
the initial destabilizing energy. Even widi tliis criticism, 
the above seiniempirical approach allows deduction of 
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FIGURE 7.11 Isudensitrace of portian of the morning trail of 
Oetobcr 16, 1969, 50 sec after release, showing the division into 
laminar and turbulent layers, Altitudes are shown in kilometers. 


many reasonable properties of the atmosphere. In particu- 
lar, ^e model counters the objections raised to the exis- 
tence of turbulence in what ostensibly is a highly stable 
region; the presence of turbulence itself tends to de- 
stabilization by modification of the temperature proBle, 
Up to this point, major emphasis has been placed on the 
fundamental parameter e, the rate of dissipation of turbu- 
lent energy. There is an equally important although not so 
easily deftned parameter, K„ die coeflicient of turbulent 
eddy diflusion in the vertical, which is the transport 
parameter incorporated in all one-dimensional models of 
the chemistry and constituents of the lower thermo- 
sphere. For some time, there was considerable discrep- 



FICURE 7.12 Theoretical morning temperature structure de- 
duced from the profile of Figure 7.11. 

oncy between the difltision coefficients calculated from 
the growth of rocket-released tracers and those inferred 
from measurements of diffusive separation and atmo- 
spheric heat-budget calculations. Widi the discovery that 
the turbulence responsible for the enhanced diffusivity of 
chemical releases in the lower thermosphere was highly 
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FIGURE 7.13 An esample of temperature structure in the 
lower stratosphere (Mitchell and Prophet, 1969). 
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FIGURE 7.14 Thf turbulent vetlkol cddV'dilTusirTD I'liulRclent 
priifilri deduct'd frnm the stiucturr on the October 1969 re- 
leases. d indicates the coefRcient deduced from simultaneous 
radio-meteor observations at Adelaide (35* S), apprusimatelydSO 
km southeast of Woumera. 


anisotropic with horizontal scales ten times those of the 
vertical, this discrepancy was readily explained. Since the 
time constant t* for the observed onset of turbulence as 
used by Rees ef al. (1972) is characteristic of the small- 
scale, isotropic turbulence spectrum, Lloyd e> ol, (1972) 
used the c values thus determined to estimate the vertical 
diffusion coeflicient. By an extension of their temperature 
profile mudifleatiun model, they equated the vertical de- 
stabilizing influence of the turbulence to the stabilizing 
influence of the mean-temperature profile to obtain a 
vertical diffusion coeflicient 



where g is the acceleration due to gravity, To is the 
temperature atallihidez, dro/dz if the undisturbed mean 
temperature gradient, and P is tlie adiabatic temperature 
lapse rate, 9 K/km. The above relationship is based on the 
fact that as lung as the vertical temperature gradient 
remains greater than the adiabatic, there is an inhibition 
of vertical transport by vertical motion. It is interesting to 
compare this relationship with that determined independ- 
ently by Lilly et al. (1973) for tlie lower stratosphere. 
The constant fi above was dclcniiincd from turbulence 
theory to be 10. Lilly arrived at the value of 1/3! This 
discrepancy is not so serious as it appears al first sight, 
since Lilly’s temperatures and temperature gradients 
were the measured values, leading to a considerably 
smaller denominator than that produced by the use of the 
undisturbed mean temperature and gradient values by 


Lloyd el al. The variation of K, with height for the 
October 1969 releases is shown in Figure 7.14. Tlie 
theoretically deduced upper limit for eddy diffusion, 
based on tlie heat flux model of Johnson and Gottlieb 
(1970), is shown for comparison. Note tliat the maximum 
value of eddy diffusivity estimated by Johnson and 
Gottlieb is the integrated global maximum averaged over 
all seasons and can be exceeded by a particular measure- 
ment if there is any temporal, latitudinal, diurnal, or 
seasonal variation. Note also tliat all tlie turbulent inten- 
sity and diffusivity proflles presented here reflect the 
apparent sharp cutoff in turbulence as observed on chem- 
ical trails. Witli the time constant in onset of turbulence 
increasing so rapidly at and above the turbopaiise, most 
trails ate not observable long enough for breakup to 
occur. However, ubsen'ed diffusivities between tlie tur- 
bopause and approximately 130 km ore usually greater 
than molecular, in agreement With the Johnson and 
Gottlieb calculations. 

The only two sets of data so far produced that are 
amenable to analysis in terms of month-by-month varia- 
tion of * or K, were determined from meteor trail shear 
measurements made at Adelaide (35° S) by Roper (1966a) 
and McAvaney (1970). The variation of the monthly mean 
K, at 93 km is shown in Figure 7.15. Also plotted in this 
figure are the results deduced from a far more rigorous 
treatment of the 1961 c data by Zimmerman; some mod- 
ifications of absolute values is evident, but the overall 
variation, with equinoctial maxima, remains. There is a 
real difference between the 1961 and tlie 1969 data that is 



FIGURE 7.15 The southern hemisphere miiilalilude variation 
of vertical eddy diffusivity deduced from two years of radio, 
meteor wind-shear observations. Zimmerman's values result 
from a more rigorous analysis of the same 1961 data. 
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ni)l one of loraJion, interprobilion, or lochnique. TTie 
hiphcr values im'usurvd in IU6H may be a eonsequenee of 
the higher solar oetivily Ihul year. 


7,4 SUl'RCES OF TURBULENCE ENEACY 

The fuel that random internal gravity waves are the source 
of the turbulenee energy in the lower thermosphere has 
already been mentioned, tVliile this subject is covered in 
nmre detail in Chapter 6, one further correlation is perti- 
nent to this discussion. In looking for possible sources of 
turbulence in the meteor region, comparisons were made 
with the magnitudes of the prevailing and tidal winds and 
shears as meosured simultaneously by the radio-meteor 
melliod. miile no apparent relationship existed between 
the turbulence parameters and the prevailing and 
semidiurnal winds, the seasonal variation ofttie amplitude 
of the diurnal oscillation was highly correlated with the 
turinileni intensity. The monthly means of the amplitude 
of tl>e diunial Ude for several years of observatian are 
shown in Figure 7.16. Note that equinoctial maxima in 
the diurnal tidal wind ainplitudc are not observed in all 
years and, therefore, that the measured variation in turbu- 
lent intensity may not be characteristic of all years. 

The global-scale diunial tidal wind does not produce 
turbulence directly but by a cascade process in which the 
tidal wind becomes unstable and generates in situ a 
spectrum of random internal gravity waves. This instabil- 
ity in the diurnal tidal wind may come about either 
because its amplitude becomes large or because of non- 
linear interactions with gravity waves generated below 
and propagating upward through the lower thermo- 



FICURE 7.16 The seasonal variation of the southern hemi- 
sphere midlutihide diurnal tidal wind energy per unit mass de- 
termined from radio-meteor data. Equinoctial maxima are not a 
feature of all years. 


sphere. Such gravity tvaves propagating from below can 
themselves be a direct source of turbulent energy. Tlius 
the dominant feature of turbulence in the lower thermce 
sphere, even if it is present at all times, will be the 
intermittency of its intensity. Little is known of the role 
played by large-scale motions such as planetary waves in 
the stability of tliis region of tlie atmosphere. In fact, since 
practically all measurements of turbulent intensity have 
been made at middle latitudes, the properties of tlie 
global turbopause and tlie influence of the turbopause on, 
for example, the hemispherical asymmetries in minor 
constituents measured globally at satellite altitudes can 
hardly be estimated at this juncture. Our knowledge of 
the temperature structure of the turbopause region at 
scales less than 1 km. which is crucial to the understand- 
ing of turbulence ond diffusivity, is woefully inadequate. 
Present measurement techniques are quite incapable of 
producing such detail at these altitudes. 

In addition to the modirications to the temperature 
prollle already discussed, the dissipation of turbulent 
energy produces heating of tlie ambient atmosphere. Fur 
the mixed atmosphere of mean molecular weight 29, a 
rate of turbulent dissipation ( of 1 W/kg produces heating 
at a rate of 85 K per day. The top scale of Figure 7.6 gives 
the heating rates appropriate to the inferred turbulent 
dissipation rates. In tbe light of the considerable variabil- 
,iy indicated by the measurements presented here, even 
more emphasis must be given to the rate of dynamical 
heating of the lower thermosphere, as proposed by Hines 
(1965). The dissipation of wind energy at these altitudes 
may, at times, give rise to local heating rates in excess of 
those due to the solar input, which has usually been 
considered to be the major source of heating in this 
region. The relative importance of turbulence in mixing 
and dynamical heating has been summarized by Johnson 
(1975). 


7.5 CONCLUSIONS 

Uhile the turbulence in the lower thermosphere is iso- 
tropic to scales of a few hundred meters, the transition 
from turbulence to the nondissipative scales of the 
gravity-wave spectrum is gradual, and therefore horizon- 
tal diiliisivity is greater than vertical diffusivity. Based on 
diffusivity criteria alotie, one would say that the turbu- 
lence is anisotropic, with a vertical scale of the order of 1 
km and a horizontal scale of a few kilometers. The turbu- 
lent intensity is intermittent in both space and time, with 
large diurnal and significant seasonal and possible solar- 
cycle variations. For this reason alone, it is essentia! that 
simultaneously measured in situ values of atmospheric 
parameters be used in any attempt at meaningful com- 
parisons. Since practically all measurements of tiir- 
bopause altitude and turbulent intensity have been made 
at middle latitudes, and those low and high latitude 
measurements that have been made have not been coor- 
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diimli'd with fiinuiltani'ims obM-rvations (‘Istwhi'n.*, varia- 
(iniis with latitiidf art- latei'l>’ unknown. However, il has 
been established by roekvl-yrenade and falling'Sphere 
measurements that the winter polar mesopause is warmer 
diati that at the summer pule and that in the absenee of 
solar input a signifleant atmcspherie dynamical heating 
iouree is responsible. As has been shown, the turbulent 
dissipation of wind energy in the 90- to 125>km height 
range is signilieant and should therefore be included in 
any model of tlie tliermuspheric heat budget. 

\Miile knowledge of the veitical diffusivit)' in the lower 
thermosphere is vital to the understanding and modeling 
of thermospheric constituenis, the rule played by the 
global turbupause in the dynamics of the thermosphere 
has yet to be determined. High-resolution photography, 
with a frame rate of at least one every 2 sec, is able to 
produce data on turbulent intensities and vertical dif- 
fusivities from riicket vapor trails, but such measurements 
are highly localized in time and space. A recent develop- 
ment, the use of a high-flying aircraft as a camera plat- 
form, overcomes two ground-based camera problemsi at- 
muspberichaze and clouds are avoided, and photographs 
can be taken in the daytime. Daytime releases of lithium 
can also be observed from the ground using narrow band 
filters and electronic scanning systems (recording on 
video tape), but interpretation of the dispersion of the 
highly energetic release in terms of turbulence parame- 
ters is difficult. Vertical dilTusivities and turbopause 
altitudes can be inferred from rocketbome mass-spectrom- 
eter measurements, regardless of hour of day or 
cloud cover but with location limitations similar to the 
vapor-trail method (both require reasonable range 
facilities). Multistation meteor wind radars can provide a 
continuous measurement of the turbulent intensity just 
below the turbopau-se, but only a few stations, all in 
midlatitudes, are currently capable of this type of opera- 
tion, and none operates continuously. Incoherent-scatter 
radars are able to measure tlie temperature structure in 
the neighborhood of the turbopause with about 2-km 
height resolution, but these measurements are subject to 
constraints similar to, but even more severe than, those of 
the radio-meteor technique. A suggestion by Beneze 
(1970) that an ionosonde can be used to measure tur- 
bopause altitude warrants further investigation, since a 
global network of ionosonde stations is already in exis- 
tence. 

A proper understanding of the nature of the turbulence 
in the lower thermosphere ntquires a knowledge of the 
temperature profile with better than 1-km (preferably 
100-m) resolution — a resolution that is not technically 
feasible at tliis time. Nevertheless, global variations in 
turbupause altitude and intensity can be determined 
using well-established techniques but only with interna- 
tional cooperation. Simultaneous global observations 
could be coordinated through the Middle Atmosphere 
Program (MAP) of the Special Committee on Solar Terres- 
trial Physics (SCOSTEP) and the International Meso- 


spheric and Ionospheric Structure Parameter Interaction 
program (siispi) proposed by the Soviet Union. Particular 
emphasis is being placed by these programs on the en- 
courogement of the Arctic, equatorial, southern hemi- 
sphere, and Antarctic observations so badly needed for 
global coverage. 
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CHAPTER IV 

Mesospheric and Lower Thermospheric Dynamics 

by 

R. G. Roper 

School of Geophysical Sciences 
Georgia Institute of Technology 
Atlanta, Georgia 30332 U.S.A. 

INTRODUCTION 

The region between 60 and 85 km altitude (the mesosphere) has often 
been referred to as the "Ignorosphere," In that It Is above the heights 
sounded by conventional meteorological rockets, and below the E region, 
which has been monitored for decades using radio techniques. 

The advent of the rocket released vapor trail in the late 1950’s did 
little to increase our knowledge of the mesosphere, in that these sodium 
vapor and trimethyl aluminum releases were concentrated In the 85 to 240 km 
altitude range, with the exception of a few titanium tetrachloride releases 
in the mesosphere. 

Since the early 1950's, a considerable amount of data on neutral atmo- 
sphere motions In the height range 80 to 100 km has come from the meteor 
radar technique. However, with the notable exception of those stations 
operated by the Hydrometeorological Service of the U.S.S.R. (which, unfor- 
tunately do not routinely measure height, but consider all echoes as coming 
from a mean echo height near 95 km), the continuous operation needed to 
Infer the synoptic meteorology of the region has been rare. The only attempt 
at such operation producing height/time profiles of winds between 80 and 
100 km has been carred out at Atlanta (34®N, 84®W), and winds from August 
1974 through December 1977 have been published by Roper (1978). While 
these results are valuable, being from a single station they provide only 


a crude insight into the synoptic meteorology of the meteor region. 

PREVAILING WINDS 

While considerable variability exists in the meteir winds measured by 
various groups, global circulation systems with seasonal variability, re- 
producing year after year, have been found. Analysis of several years 
northern hemisphere radio meteor wind data by Minina, et ai. , 1977, and 
Illichev and Portnyagin, 1977, shows that the pressure field at meteor 
heights in winter is similar to the baric profile of the underlying atmo- 
sphere. During the summer, however, the lower thermosphere differs consid- 
erably from the baric profile of the mesosphere, stratosphere, and tropo- 
sphere. In winter, at high and mid latitudes a cyclonic vortex with its 
center near the pole is observed; in summer, at latitudes greater than 65®, 
anticyclonic motion prevails. Cyclonic motions prevail through the year 
at mid-latitudes, with anticyclonic circulations in the subtropics. On the 
average, the spring reversal of the circulation in the meteor zone occurs 
before that in the stratosphere--thus continuous monitoring of mesopause 
circulation can be used as a predictive tool in inferring the circulation 
of the upper stratosphere. 

To date, the meridional circulation, which is of comparable strength 
to the zonal at meteor heights, has received little detailed attention. 
However, there does seem to exist a global system of motion from summer to 
winter pole, which agrees with the warm winter/cold summer mesopause inferred 
from high latitude rocket grenade temperature soundings. 

One outstanding feature of the prevailing wind motions as deduced at 
mid-latitudes (see, for example, El ford, 1974) is the fact that the 85 km 


(mesopause) level Is an obvious boundary between the mesospheric circulation 
below, and the thermospheric circulation above. 

MEAN VERTICAL MOTIONS 

One of the most puzzling features of winds measured In the 80 to 100 
km height range Is the apparently large (meters per second) mean motion 
which Is Inferred from both Incoherent scatter and radio meteor wind obser- 
vations. While the large magnitude of these winds Is still In question, 
the sign of these motions, at least as Infisrred from the quasi -continuous 
observations over Atlanta, Is consistent with the general circulation 
(Dolas, 1979). 

TIDAL WINDS 

The most outstanding feature of the motions at meteor height Is the 
prevalence of solar tidal oscillations. On the average, more wind energy 
resides in tidal motions than in the prevailing circulation in the lower 
thermosphere. 

While the Irregularity of the diurnal tide, particularly in the north- 
ern hemisphere, has hampered the investigation of this tide as a global 
phenomenon, recent work by Mathews (1976) using the Arecibo Observatory in- 
coherent scatter facility, and by the French with their portable meteor wind 
radar at Ramey, Puerto Rico (results as yet unpublished) has shown the pres- 
ence of a stable, large amplitude (~50 m/sec) diurnal tide at the equinoxes. 

In the southern hemisphere (Adelaide, 35®S, 139°E), the diurnal tide 
has a larger average amplitude than is found at northern hemisphere mid- 
latitudes, but both phase and amplitude are highly variable In both hemi- 
spheres. 


The semidiurnal tide, by comparison, exhibits relatively constant 
phase, but, again, great variability in amplitude. This variability can 
be explained by variations in the stratospheric/mesospheric wind fields 
through which the tidal energy is propagating to reach the mesopause. The 
variability has time scales related to the planetary wave time scales in 
the lower atmosphere, and further understanding of these variations can only 
come from long period (p«jrf€!rably continuous) observation by a network of 
stations widespread in both latitude and longitude. 

STRATOSPHERIC WARMINGS 

It has been known for some time that the sudden warmings of the polar 
stratosphere which occur during some northern hemisphere winters affect the 
circulation at mesopause levels at high latitudes. More recent measurements 
at mid-latitudes (Atlanta, 34“N, 84®W) have shown a direct cause/effect re- 
lationship, with major stratwarms causing zonal wind reversals. There are 
indications from the limited amount of data available from the southern hemi- 
sphere that northern hemisphere winter polar stratwarms may be a global 
phenomenon at the mesopause level. 

Stratospheric warmings also affect tidal amplitudes and phases, but 
with a lack of cos^tlnuous height/time profiles available globally, reports 
of such effects have been for the most part descriptive of particular loca- 
tions only, and somewhat confusing in their interpretation. 

MODELS 

Two significant semiempirical models of stratospheric/mesospheric/ 
lower thermospheric winds have emerged in the last decade. Groves (1970) 
took data from 1000 rocket launches with ejected sensors, 127 rocket gre- 


nade experiments, and 230 experiments with rocket released vapor trails and 
clouds, and constructed monthly mean zonal and meridional circulation pat- 
terns. This model, together with monthly mean tropospheric data, upper air 
global temperature maps, radio meteor winds, and the Jacchia thermospheric 
model have been Incorporated In the Global Reference Atmospheric Model of 
Justus, et ai, 1976. This empirical model generates latitude, longitude, 
and altitude dependent monthly mean values for pressure, density, tempera- 
ture and winds from surface to orbital altitudes. 

Theoretical models of atmospheric tides have recently been reviewed by 
Forbes and Garrett (1979). While there Is a measure of agreement between 
the latest theoretical models and observational data, since the tides are 
a global phenomenon, questions as to standing versus travelling waves, and 
latitudinal variation will not be resolved until a global network of synop- 
tic measurement capability Is established. 

An additional problem In resolving theory and observation centers 
around the variability encountered in all meteorological phenomena. The 
presence of other motions besides those of primary interest, in particular, 
the aliasing of tidal observational data by long period (planetary) waves, 
and the interaction between the tides and the shorter period random internal 
atmospheric gravity waves, makes a widespread, continuously operating moni- 
toring network a necessity. An ideal radar for deployment is that described 
by Aso, e;t al_* (1979). 

COMPATABILITY OF MEASUREMENT TECHNIQUES 

While early attempts to compare winds measured by different techniques 
(e.g., rocket released vapor trails, radio meteor winds, airglow drifts. 


and partial reflections radio observations) were singularly unsuccessful, 
more recent comparisons, in which allowances were made for sampling differ- 
ences, in particular integration Intervals in both height and time, have 
resolved these inconsistencies (see, for example, Felgate, et 1975; 
Geller, e^al_. , 1976; Vincent, et a],. , 1977; and Hernandez and Roper, 1979), 
As yet unpublished comparisons between the French meteor radar located at 
Ramey, Puerto Rico, and Arecibo incoherent scatter facility, have shown 
excellent agreement. One outstanding achievement at Arecibo has been the 
extension of the lower limit of reliable wind determination by the incoher- 
ent scatter radar down to as low as 55 km. 

The radio meteor wind technique is capable of continuously monitoring 
the dynamics of the 80 to 100 km altitude range, with a lower limit of ob- 
servable variation of a few hours, and a height resolution of some 2 km. 
Specialized highpowered radars, such as that operated by the University of 
Illinois, can monitor variability of time constant less than one hour, but 
the logistics of highpowered operation preclude continuous sampling beyond 
intervals of a few days. The diurnal variability of meteor influx renders 
the error of wind determination greater at dusk than at dawn, but most facil 
ities achieve a sufficiently high echo rate at dusk to make these results 
still statistically significant. 

The partial reflections drift experiment is capable of sampling, with 
a time resolution of minutes, and a height resolution of 3 km, a height 
range from 60 kilometers to 110 km in daylight hours, but returns are reli- 
ably received only from 90 km up at night. The incoherent scatter technique 
suffers similarly from the diurnal variation of mesospheric ionization. 

Rocket vapor techniques, while yielding excellent instantaneous snap- 
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shots of the stratosphere/inesosphere/ thermosphere (height range governed by 
type of vapor released, and time of day) are prohibitively expensive if 
synoptic data is required. Their dependence, together with airglow observa- 
tions, on cloud cover for observation is also a significant drawback. 

Winds can be inferred in the stratosphere, mesosphere, and lower thermo- 
sphere by application of the thermal wind equation to satellite temperature 
soundings. Unrortunately, the large tidal winds at upper mesospheric heights 
and above, preclude the use of quasi-geostrophic approximations— the inference 
of viable winds from satellite temperatures at these altitudes is highly 
doubtful . 

A further groundbased technique worthy of mention is the VHP radar 
technique (Rottger, et §2.., 1978). While not yet in routine operation, 
these VHP radars have demonstrated the capability of measuring winds and 
turbulence in the troposphere, stratosphere, and upper mesosphere. 

Also, the VLP observations of Schminder and Kurschner (1979) warrant 
further investigation, since they appear to provide continuous, reliable 
observations of mesopause dynamics. 

AVAILABILITY OP DATA 

Beginning with the formation of GRMWSPM, the Global Radio Meteor Wind 
Studies Project of lAGA, in 1970, followed by the URSI/IAGA Cooperative 
Tidal Observations Program (CTOP), which added incoherent scatter radar ob- 
servations, attempts have been made to coordinate observations (see Roper 
(1979)). Some of these results were published in the August, 1978, issue 
of the "Journal of Atmospheric and Terrestrial Physics." 

While data is available from individual observers, as yet formats and 
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analysis techniques have not been standardized, and no central depository 
enlsts. The Soviet Union Is moving toward having all their radio meteor 
wind data available through the Uorld Data Center In Moscow. The mechanics 
of such archiving, both In the Soviet Union and elsewhere, are to be dis- 
cussed at a meeting of IA6A Division V, Working Group 2 (Meteor Observator- 
ies) at the lUGG Assembly In Canberra, Australia, In December 1979. 
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DYNAMICS OF THE EQUATORIAL MESOPAUSE 
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ABSTRACT 

Until recently, our view of the neutral atmosphere dynamics of the equatorial mesopause 
has been restricted to "snapshots" resulting from either one time campaigns, or Intermittent 
observations. Rockets flown from Brazil, the "big gun" launches from Barbados, meteor wind 
radar data from Somalia and Jamaica, and VHF radar observations from Peru, while producing 
much information on the nature of motions with periods of up to a day, provided only meager 
details of the synoptic seasonal meteorology of the region. More recent observations have 
resulted from the extension to lower altitudes of the Incoherent scatter technique at Arecibo, 
the Installation of a partial reflection drift station at Townsville, Australia, and the 
establishment by the French of a meteor wind radar at Ramey, Puerto Rico. The Incoherent 
scatter and partial reflection facility operate In the "campaign" mode while the French 
meteor radar Is operating continuously. 

Because of the limited amount of synoptic data available, this review concentrates on 
results produced by those techniques most likely to be used for synoptic observations In the 
future. However, since even these promising techniques may not be able to be used In other 
than a campaign mode (for reasons ultimately financial), a rationale for their use In winter 
months (on the basis that northern hemisphere winter polar stratospheric warmings are probably 
a global phenomenon at mesopause altitudes) Is presented. 


INTRODUCTION 

Many observational techniques have been 
used to measure winds at mesopause heights 
over the tropics - chemical trails from pro- 
jectiles launched by the "big gun" on 
Barbados (Murphy, et al., 1966), rocket 
grenades over Natal, Brazil (Groves, 1974), 
radio meteor measurements over Somalia 
(Babajamov, et al., 1970), and Jamaica 
(Alleyne, et al., 1974, Scholefield and 
Alleyne, 1975), VHF radar observations 
(Woodman and Guillen, 1974), and Incoherent 
scatter observations (Mathews, 1976; Harper, 
et al., 1980). More recently, a partial 
reflection drift station has been set up at 
Townsville, Australia by the Department of 
Physics of the University of Adelaide 
(Dr. R. A. Vincent), and a meteor radar by 
the French National Center for Telecommuni- 
cations Studies (Dr. M. Glass) at Ramey, 
Puerto Rico. The various techniques have 
been reviewed most recently by Woodman 
(1977) and Evans (1978). 

Woodman places particular emphasis on 
equatorial measurements and the need for 
more observational data. His paper at this 
meeting concentrates on small scale struc- 
ture (gravity waves and turbulence), and 
so this present work will confine itself to 
tidal and planetary waves and the prevailing 
winds. 

OBSERVATIONS 

To date, all observations of winds at 


the mesopause level over the equator have been 
of the "snapshot" variety, with a few days of 
observation, very occasionally repeated sea- 
sonally. (The same can also be said for other 
latitudes, with the exception of the meteor 
radar network In the U.S.S.R. operated by the 
Hydrometeorological Service, which operates 
continuously, but unfortunately does not mea- 
sure wind structure with height.) These 
"snapshots" have revealed considerable vari- 
abllitjf In the winds, but with Indications of 
predominant easterlies at the equator. In 
contrast to the predominantly westerly flow 
observed at mid latitudes. 

RADIO METEOR WINDS 

An example of one of these observations, 
obtained by a collaborative effort between 
the French National Center for Telecommunlca- • 
tions Studies, and the Georgia Institute of 
Technology, using the French meteor radar 
located at Ramey, Puerto Rico (18°N), Is pre- 
sented In Figure 1. The zonal wind component 
of Figure la Is predominantly easterly (wind 
vector directed toward the west), for the 
first four days of observation, but a strong 
westerly intrusion then develops from above. 
The Intrusion weakens and ascends toward the 
end of the observational period. While this 
phenomenon can be associated with the simul- 
taneous spring reversal of the zonal strat- 
spherlc circulation. It may also be evidence 
for the mean flow acceleration by tidal 
interaction proposed by Miyahara (1979) - as 
can be seen from Figures lb and c, both the 
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Figure 1. Radio meteor winds from 85 to 105 
km for the period March 17-24, 1979, over 
Puerto Rico, a) zonal wind; b) 24 hour 
tidal amplitude; c) 12 hour tidal amplitude. 

diurnal and semidiurnal tidal amplitudes 
have large amplitudes coinciding with the 
reversal . 

One phenomenon which has received 
considerable attention at middle latitudes 
is the 2 day wave, which has been observed 
in week to 10 day meteor radar winds in both 
the northern and southern hemisphere. By 
using a low pass filter to remove periods 
shorter than 30 hours from 8 days of meteor 
wind radar data obtained at Ramey in August- 
September 1977, Glass (private communica- 
tion, 1978) produced the results shown in 
Figure 2. The presence of a wave with a 
quasi -2 day period is well illustrated. 

INCOHERENT SCATTER WINDS 

Some recent results from Harper, et a1. 
(1980), taken using the Arecibo incoherent 
scatter radar from September 1 to 14, 1977, 
are shown in Figure 3. The 14 noontime 
orofiles illustrate the variability encoun- 


LONB PERIOD VflRIfiTlONS OF THE ZONAL MUa 
RANEY -flU0-6EPT 1977 



ALTITUDE 

imi) 

tor 


83 


85 


82 

89 


86 


9 1977 


Figure 2. Radio meteor winds - long period 
variations. 
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Figure 3. Meridional winds from the Arecibo 
incoherent scatter radar (from Harper, et a1., 
1980). 


tered in winds at the mesopause level. Note 
that these are meridional winds and that 
speeds of tlOO m/sec occur - zonal and merid- 
ional wind speeds are comparable at altitudes 
above the mesopause, in contrast to the pre- 
dominantly zonal flow in the stratosphere and 
mesosphere below. 

METEOR/ INCOHERENT SCATTER WIND COMPARISONS 

Mathews, et al. (1980) presents several 
comparisons between meteor winds from the 
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Figure 4. Location and sampling regions of Ramey meteor and Arecibo incoherent scatter radars 
(from Mathews, et a1., 1980). 
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Figure 5. Comparison between Arecibo incoherent scatter radar winds (dashed lines) and Ramey 
meteor radar winds (full lines), August 3, 1978 (from Mathews, et al., 1980). 

French radar at Ramey and winds from the using this technique. However, there are no 

incoherent scatter radar at Arecibo. The incoherent scatter winds below 100 km at 

locations and "zones of echo acceptance" of night, whereas the radio meteor technique 

the two instruments are shown in Figure 4. produces winds over the 80 to 100 km height 

One comparison, from 0900 to 1600 hours on range both day and night. 

August 3, 1978 (Figure 5), shows quite good 
agreement between the two techniques, with 

the incoherent scatter winds showing more PARTIAL REFLECTION DRIFTS 

structure because of better height resolu- 
tion. Between 1000 and 1400 hours, the inco- This relatively inexpensive technique is 

herent scatter winds are significant down to capable of measuring winds between 60 and 100 

60 km - a real breakthrough when one remem- km by day, and 90 to 100 km by night. The 

bers that just a few years ago 110 km was first partial reflection drifts experiment 

the lower limit of reliable wind measurement operating on a routine basis was installed at 






Townsville, Australia (20°S) by Dr. R. A. 

1 Vincent of the University of Adelaide, 

I Australia. This equipment Is capable of con- 

I tinuous operation, but for financial reasons 

I Is run In the campaign mode. Results from 

{ this experiment are eagerly awaited. 

I VHP SCATTER RADARS 

I The most recently developed tool cap- 

I able of continuous monitoring of mesopause 

I dynamics Is the VHP Scatter Radar. These 

I radars have the capability, shared by the 

I partial reflection drift technique, of mea- 

i suring periodicities In wind profiles down to 

I the Brunt Valssala, and therefore represent 

I powerful tools for gravity wave and turbu- 

I lence studies. In addition to planetary waves, 

I mean winds and tides. Details of this tech- 

nique are to be presented at this meeting by 
Ben Baisley. 

I AIRGLOW DRIPTS AND TEMPERATURES 

[ While unable to provide continuous mea- 

surements (nighttime only, no clouds), alr- 
i glow measurements can supplement our know- 

V ledge of winds (see, for example, Hernandez 

and Roper, 1979) and provide temperature data 
of use In modeling mesopause dynamics. 

POLAR STRATWARM EPPECTS AT MESOPAUSE HEIGHTS 

The necessity for continuous monitoring 
of mesopause dynamics Is illustrated In 
I Figure 6 (from Roper, 1978). Without 



Figure 6. Comparison between zonal ly aver- 
aged satellite radiance data and the zonal 
wind at meteor heights over Atlanta (from 
Roper, 1978). 


continuous data, the association between the 
polar winter stratospheric warming and the 
reversal In the zonal wind at mesopause, alti- 
tudes would not be nearly as convincing. 

While this data Is for mid-latitudes (Atlanta, 
34°N), the effect Is so pronounced that It 
would not be surprising If It extends as far 
south as Puerto Rico (18°N). 

CONCLUSIONS 

No mention has been made here of 
satellite remote sensing to monitor meospause 
dynamics. Because winds at these altitudes 
are ageostrophic, conversion of temperature 
fields to winds Is well nigh Impossible. For 
many years to come, and certainly through the 
Middle Atmosphere Program period, 1982-85, 
ground based radio techniques will be our 
primary source of data on mesopause dynamics 
at all latitudes. 
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Abatract 

Hagnitudes arid horiapntal vercicf.1 scales 
of gravity waves have been neasared and uied in a 
formula developed by tiinos to ettimate resulting 
eddy diffusion coefficients. Values between 90 and 
110 km agree well with those determined from energy 
flux measurements from chemical cloud releases. 
Values between 23 and AS km altitude are in agree* 
ment with high eddy diffusion coefficient! in the 
upper stratosphere recently determined by Uofsy. 
Comparison with other eddy diffusion estimates is 
also given* 

1. Introduction 

A recent study of upper atmospheric data In 
the region above 25 km has been conducted for the 
purpose of resolving small scale variations which 
may be associated with internal atmospheric gravity 
waves. The method used* known as the daily differ* 
ence' technique, was developed for resolving small 
scale irregular variations from limited amounts of 
data (i.e. of insufficient number for direct 
Fourier analysis) ^ The form of the method used in 
this study is a slight extension of the method orlg* 
inally developed by Woodrum and Justus^^*^). The 
magnitudes and horizontal and vertical scales of 



Fig. 1. Gravity wave wind magnitudes 


the small scale perturbations, assumed to be pro* 
duced by gravity waves, have been measured by this 
technique. The detsils of the daily difference 
metnod and the results of these magnitude and wave- 
length measurements hove been presented by Justus 
and Woodrum^^^ * 

Figure 1 shows the computed height variation of 
the gravity wave winds averoged over all seasons and 
locations from Meteorological Rocket Network, fall* 
Ing sphere, meteor radar, grenade release, and chem- 
ical release data. Horizontal and vertical wave* 
length results are shown in Figure 2. Data points 

Vertical Wavelength, km 



Fig. 2. Gravity wave horisontal 
and vertical wavelengths 

in Figure 2 are plotted at the midpoints of the 
altitude intervals over which they were evaluated. 
These averaging Intervals are indicated on the 
right hand side of the figure. Figure 2 shows that 
horizontal wavelengths are approximately 10 times 
larger than vertical wavelengths throughout the 
altitude interval measured. 

2. Eddy Diffusion Coefficients 

has developed a relation for the 
eddy diffusion coefficient Kp corresponding to 
"offset condition" dissipation of internal atmos- 
pheric gravity waves. The offset condition is the 
situation in which viscosity and thermal conduction 
Just offset the normal (i.e. no*visco8lty) 
exponential increase of gravity wave amplitudes, 
yielding constant amplitude. The formula. 
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cxpratMd In ectw of horliontal and yarctdal wave- 
langtha >1 and preaaura acale height Hi and 
Brunt-VataaU period tg, la glvan by 

Kp • O.OW 1^* <t, 8)"^ «) 

Ihe waveiength valuea froa figure 2, together with 
the 1962 Uf S. Standard Ataotphere valuea for t 
and H have been uaed to evaluate Kp froa equation 
(!)■ The raeulta are ahown In figure 3. The aoltd 
aquare data polnta correepond to the wavelength 
average value data polnta of Figure 2. The aolld 
curve repreeenta valuer coaputed froa the Interpo* 
latlon curves of wavelength valuea In Figure 2. 

Spaa eatrapolaelon of wavelength valuea waa 
eaployad, but not beyond the altitude llaite of the 
original data from which the wavelengtha Were 
detemlned. The portion of the aoltd curve between 
23 and 3$ ka in Figure 3 le froa the daehod eetra- 
polatlon curve in Figure 2, end the dashed portion 
of the present results curve In Figure 3 Is froa 
the dotted estrapoletlon curve (which goes to eero 
at the Earth's surface) In Figure 2, In the case 
of << 1^1 equation (1) reduces to 

Kp • 0,014 (tg H Xg)*‘ (2) 


Thus I since Xg « ^0 values vary approxi* 

ttateiy as ths -fourth power or X^p but inversely with 
Xgi It is for this reason that snoll changes in Xg 
in the 25 to 35 ko altitude interval (where X^ is 
Bfflall and stnall changes represent large percentage 
changes) make considerable difference in the values 
of Kp computed from equation (1). 



Figs 4s Gravity wave kinetic energy density 

Eg exp<»e/7'e8) I These vslues of "kinetic energy 
density scale height" which are not significantly 
larger than attnospheric density scale height* indi- 
cate; thot the tctuel rote of gravity wave dissipat'^ 
ion Bore nearlv corresponds to the required oftset 
condition than the pV^ • constant situation without 
viscous effects • 

3s Comparison With Other Results 



The gravity wave magnitude results of Figure 1 
show that the gravity wave amplitudes arc* on the 
average, Increasing with height, contrary to the 
assumed offset conditions necessary for application 
of equation (l)i However, as shown in Figure 4, 
there is considerable lose of gravity wave kinetic 
energy density with heights Between about 25 and 
60 km the kinetic energy density E ■ pV^/2 varies 
as £ • Eq cxp(-z/9.4) with a constant and z in 
km. Between about 70 and 130 km E varies as 


by Eg • Kp \ curve labeled 
resuitfi ft/o (heir computat- 


Flgure 3 also shows values and curves of eddy 
diffusion coefficients which have been measured or 
used below 110 km. The dashed curve denoted Webb'^' 
was also based on earlier work of Lettau^'^ in the 
neutral atmosphere altitudes. The curve labeled 
Justus'®' was determined from measured rates of tur- 
bulent k'l netic energy production Cg by wind shear* 
and BU!ar wind ahear S, with Cg assumed to be re- 
lated to Kn ami S bv f ■ ■ 

Johnson ano Wilkins' 
ion of the maximum eddy diffusion c ■ aistent with the 
amount of heat available for downward eddy 
port* The data points labeled Gudlksen et al. ^ 
were determined by these authors from their obscr- 
vacions of nuclear explosion released tungsten 165. 
The data pointn labeled Kellogg(H) were computed by 
the present author from Kellogg’s reported values of 
turbulent velocities (v^) and maximum observed rad- 
iUB (r(Qgx) umoke puffs in the 10 to 18 km height 
range. The data points are values of the product 
rgjgy, which should represent o lower limit to the 
eddy oif fusion coefficient if, as suggested by 
Kellogg, the linear growth velocity is due to in> 
Crcoslxng eddy diffusion coefficients with increas- 
ing scale. However, as discussed below* this con- 
clusion is quite likely incorrect. 


4. Interpretation of Results 


Because vertical spread within the atmosphere 
is ultimately rcstricced compared to the unlimited 
hori^'Kontal spread which Is possible* the large 
scale eddy dlflusion process is not isotropic, and 
vertical eddy diffusion coefficients are generally 
smaller than h(>rizontal eddy diffusion coefficients 


The first point to be noted with regard to the 
eddy diffusion coefficient data in Figure 3 is chat 
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(with the probably exceptloi of the values conputed 
frota Kellogg's snoke puff d(ica) these results ere 
for vertical eddy diffusion coefficients. The 
Johnson and Wilkins, Wofsy, and Gudlkssn at el. 
date all relate theoretically and/or observot* 
tonally to vertical fluxes of diffusing constlt> 
uents or heat. As s Dessure of the degree of 
anisotropy possible, it should be noted that the 
purely horieontal (aerldionul) eddy diffusion 
coefficients determined by (iudlksen et el. are 
roughly a factor of 10° larger than their purely 
vertical eddy diffusion coel’ficients plotted in 
Figure 3. 

It ie likely that the rg,aK values from 
Kellogg’s snoke puff date ovetestlnace the vertical 
eddy diffusion coefflclonis, because apparently the 
actual diffusive spread of cloud releasee on snail 
scales proceeds faster than would be accounted for 
by the large scale vertical eddy diffusion coeffic* 
lent. This can be seen fron canpariscn of diffus- 
ion velocity-scale products (l.e. V( t„.^) from 
chemical releose clouds^'*' with the eddy diffusion 
coefficient values (labeled Justus in Figure 3) 
deternined from vertical energy flux measurements 
using chenical release data. The velocity-scale 
product values are consistently at least a factor 
of ten larger. This result Is undoubtedly caused 
by the small scole diffusive spread being controlled 
by wind sheer processes. 

The present results curve in Figrre 3 should 
most precisely be interpreted as a maximum value 
which could be produced by internal atmospheric 
gravity waves without an upward decreise in the 
anplitudc of the waves. Therefore, tiis curve 
might require a slight shift to lower values be- 
fore it would represent tht profile of actual mean 
vertical eddy diffusion coefficlencs. Neverthe- 
less, it is clear from Figure 3 that the present 
results agree with high stratospheric eddy diffus- 
ion values such as those of Wofsy and cannot agree 
with the Webb-Lettau curve which has -ttratospheric 
eddy diffusion coefficients which are lower by sev- 
eral factors of ten. 

The important conclusion for theuc results, 
with regard to the assesament of environmental im- 
pact of aerospace operations, is that the vertical 
eddy diffusion coefficient in the str.itosphere, and 
especially in the upper stratosphere, is higher 
than originally believed. From the rosults pre- 
sented here it is reasonable to attribute these high 
eddy diffusion coefficients to the mixing processes 
generated by atmospheric gravity waven as they 
propagate through the stratospheric region. 
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CHAPTER VII 

Stratospheric Warming Effects 
on Kesopause Dynamics 


a) The Effects of Polar Stratwarms on the Winds at the Mesopause Level 
in Mid-Laditudes 

(Presented at the 18th Radar Meteorology Conference, Atlanta, 6a., 
March, 1978, and updated at the AMS Conference on the Upper Atmo- 
sphere, Boston, Mass., October, 1978) 


b) Prevailing Wind in the Meteor Zone (80-100 KM) Over Atlanta, and its 
Association with Mlu-winter Stratospheric Warming 

(Submitted to the Journal of Atmospheric Science) 
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THE EFFECTS OF POLAR STRATWARMS ON THE MINDS AT THE MESOPAUSE LEVEL IN HID LATITUDES 


R. 6. Roper 

School of Aerospace Engineering 
Georgia Institute of Technology 
Atlanta, GA 


1 . INTRODUCTION 

For over a quarter of a century, winds have 
been measured over the height range 80 to 100 km 
by means of radio reflections from the Ionized 
trails of meteors. Although sophisticated high 
frequency radars (usually operating In the 20 to 
40 MHz range) have been employed, not all of the 
advantages of the technique have been realized. 

In particular, very few continuous, long term 
observations have been made. In the past, the 
major problem has been an Inability to handle 
the large amount of data generated by continuous 
operation. However, the revolution In data 
acquisition In recent years had made such opera- 
tion economically feasible. 

2. THE METHOD 

The radio meteor method of wind determina- 
tion has been described at a previous AMS Radar 
Meteorology Conference (Barnes, 1972), and the 
Georgia Tech Radio Meteor Wind Facility detailed 
by Roper (1975), The Tech system Is a continuous 
wave radar (using two continuously transmitted 
carriers 720 Hz apart to determine range) de- 
signed to measure some 1000 line of sight 
dopplers and echo positions per day. However, 
since the CW technique accepts echoes at all 
ranges, aircraft reflections are a problem, at 
times reducing the rate below 300 echoes per day. 
This, together with time lost through maintenance 
of the transmitter and receivers, occasional power 
outages at the receiving site, and the need to 
average several days data to produce meaningful 
measurements of the diurnal and semidiurnal 
tides, results in an apparent smoothing of the 
data of from a few days to two weeks. In order 
to meet the publication deadline, consideration 
is given here only to the zonal mean wind; analy- 
sis of the significance of the meridional and 
vertical mean winds, and the tidal winds, Is 
proceeding. 

3. RESULTS 

Continuous radio meteor wind observations 
comnenced over Atlanta in August, 1974. Figure 
la shows the variation with height and time of 
the zonal wind 80 to 100 km over Atlanta for the 
period November, 1974 through February, 1975. In 
interpreting the structure present, the averaging 
interval details In Figure lb should be taken 
into consideration. 

An attempt was made to correlate the ob- 
served winds with the zonal mean satellite rad- 


30332 



Figure la) The mean zonal wind over the 
height range 80 to 100 km above Atlanta for 
the period November, 1974 through February, 

1975. The zonal means to which the contours 
have been fitted are averages for the measure- 
ment Intervals depicted in b). c) is the 
zonal mean radiance deduced from satellite 
radiometer measurements, and is proportional 
to the stratospheric temperature averaged 
over the 100 to 5 mb height range at 82.5°N. 

lance data at 82.5‘’N published by Quiroz et al 
(1975), The late December/early January warming 
did not produce a wind reversal at the 10 mb 
level, therefore is not characterized as a major 
warming. 

The circulation at meteor heights during the 
fall and early winter of 1974, being predominantly 
easterly (as plotted, a negative wind is a wind 
vector directed toward the west, I.e. an easterly). 
Is unusual, at least when compared to subsequent 
years (see Figures 2 and 3). About the best one 
can say from this preliminary comparison of the 
meteor winds with the single zonal mean radiance 
data curve is that there are wind reversals which 
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Figure 2. As for Figure 1 for November, 

1975 through February, 1976, except that the 
zonal mean radiance is presented for three 
latitudes, and is proportional to the 
stratospheric temperature averaged over the 
100 to 2 mb height range. 



Figure 3. As for Figure 2, for November, 1976 
through February, 1977. 


correlate with the changes in radiance. 

Proceeding to the winter of 1975-76 (Figure 
2), the zonal flow at meteor altitudes is pre- 
dominately westerly, with two minor reversals - in 
early January, and early February. Comparison 
with the zonal mean radiance satellite dat,a, this 
time available for three latitudes (Quiroz, pri- 
vate communication, 1977) shows that there are 
local maxima in the radiance curves at all three 
latitudes during these two periods. However, 
there are not any outstanding features in the rad- 
iance data - no warmings occurred curing the 
winter of 1975-76. 

The data for the winter of 1976-77 represents 
the most interesting set to date. Unfortunately, 
no meteor wind data is available fcr the first 
week in November (and the last week in October). 
However, the winds through August, September and 
mid October were consistently westerly - thus the 
change to easterly in mid NoveniUer could well be 
associated with the warming of late November. 

With the reestablishment of westerlies over the 
whole height range by the end of November, the on- 
set of the late December warming is preceded by a 
rapid waning of the westerly at the upper levels, 
with the onset of easterly flow above 90 km 
occurring a week before the winds at 80 km become 
easterly. During the last week in December, the 
satellite radiance data show a reversal of the 
latitudinal temperature gradient from that 
established by the warming. This -s accomplished 
by what one would almost call a westerly ":et" in 
the mid height range of the meteor winds. The 


subsequent warming trend in early January is not 
accompanied by a meteor wind reversal, and just so 
that we will not be tempted to think that we have 
all the answers, there is a tongue of weak eas- 
terlies which descends from 100 km to 83 km in 
late January/early February, which seems to have 
no correlation with the stratospheric temperature 
curves. 

4. CONCLUSIONS 

The continuous measurement of the wind pro- 
file over the 80 to 100 km region by means of radio 
reflections from meteor trails provides data which 
can be correlated with stratospheric temperature 
changes inferred from satellite radiance data. 

The very preliminary results presented here (later 
publications will include the consideration of the 
meridional wind, as well as the diurnal and semi- 
diurnal tidal winds) show a tendency for strato- 
spheric wanning events (which produce zonal mean 
temperature changes in the stratosphere at 40°N 
which fe only just measurable) to produce dra- 
matic changes in the circulation in the neighbor- 
hood of the mesopause at 34°N. While such changes 
have previously been expected and reported at 
higher latitudes (see, for example. Hook, 1970), 
to my knowledge this is the first report of similar 
behavior as far south as 34°N. 
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PREVAILING WIND IN THE METEOR ZONE (80-100 KM) OVER ATLANTA 
AND ITS ASSOCIATION WITH MID-WINTER STRATOSPHERIC WARMING 


ABSTRACT 

The wind data generated by an all sky, continuous wave radio 
meteor wind facility at Atlanta (34“N, 84°W) is analyzed over the 
period of August 1974 through July 1975. Zonal and meridional com- 
ponents of the prevailing wind over the height range of 80-100 km, 
at 2 km interval represent 5 to 10 day averages where the tidal com- 
ponents have been removed. 

Large southerly wind during winter and weak northerly wind 
during summer at 80-100 km altitude is consistent with other obser- 
vations and mesospheric circulation models. Evidence of wave energy 
input from the lower atmosphere has been found in the prevailing wind. 

The "minor" mid-winter stratwarm of 1974-75 is studied. Various 
phases of this warming, including a pre-warming pulse during the 
second half of November 1974, are shown to have affected the prevail- 
ing wind regime in the meteor zone, in a manner consistent with the 
latitudinal and vertical temperature compensation criteria in the 
stratosphere and above. 

The implied meriodional mean temperature gradients from the one 
year of data point to a year round cold mesopause at high latitudes-- 
the only exceptions during winter appear related to the perturbed 
stratospheric cyclonic vortex due to warming episodes. 
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PREVAILING WIND IN THE METEOR ZONE (80-100 KM) OVER ATLANTA 
AND ITS ASSOCIATION WITH MID-WINTER STRATOSPHERIC WARMING 


1. Introduction 

The general circulation in the stratosphere and lower meso- 
sphere consists primarily of ultra-long wave motions, superimposed 
on a zonally symmetric flow that is easterly (from the east) in the 
summer hemisphere and westerly in the winter hemisphere, in approx- 
imate thermal wind equilibrium with the zonal mean temperature field. 
The planetary scale waves of zonal wave numbers 1 and 2 appear to be 
vertically propagating modes that are generated in the troposphere 
by a variety of dynamic mechanisms, and then transport energy and 
momentum into the stratosphere and beyond. In addition, there is 
a semi-annual oscillation which has its maximum amplitude near the 
equatorial stratopause, and there is also a strong oscillation in 
the lower tropical stratosphere of somewhat irregular period, aver- 
aging about 26 months, called the "quasi-biennial oscillation". 
Besides these cyclic variations, there are a number of irregular 
short-term and year-to-year variations which are apparently due to 
the tropospheric forcing of vertically propagating planetary waves. 
The most prominent are the changes in the wind and temperature re- 
gimes associated with sudden warming of the high latitude winter 
stratosphere (the so-called "polar stratwarms" ) . 

On the other hand, based upon very sparsely distributed data, 
in time and location, a general profile of the seasonal variations 
of the prevailing wind in the lower thermosphere (up to 120 km) is 


becoming established. Gaov&6 (1971) finds that there are equatorial 
easterlies present throughout the year at an altitude of 95 + 15 km. 

The seasonal variation observed at low latitudes is the result of the 
slight shifting of the wind belt north and south of the equator. 

Recently tUMm et al, (1977) have analyzed the wind data from many 
years at heights of 80-100 km, primarily from meteor radar observa- 
tions made at various Northern Hemisphere stations. These data indi- 
cate large scale pressure formations of the cyclone-anticyclone tyiie, 
with similar synoptic structures as observed in the atmosphere below. 
During winter the cyclonic vortex centered near the pole dominates the 
circulation in the meteor zone at high and middle latitudes (up to at 
least 120 km; Gn.o\)QA>, 1971); in summer, the anti cyclonic circulation 
is characteristic of high latitudes (north of 65°N). The cyclonic 
circulation in the middle latitudes and the ring of sub-tropical highs 
exist year round. Thus, the winds in the meteor zone are shaped not 
only by the seasonal variations of the thermal balance at these heights, 
but also by the lower-lying layer of the mesosphere. 

The present work investigates the imprint on the general circu- 
lation at lower thermospheric heights (80-110 km) of the very large 
perturbation of the high latitude winter stratosphere which occurs 
during a polar stratwarm event. Hook (1972) using the radio meteor 
technique at College, Alaska (65°N, 140“W), and GfiZQomj and Man^on 
(1975) using the partial reflection technique over Saskatoon, Canada 
(52°N, 107°W), have studied the effects of various major stratwarms 
on the wind patterns in the 80-110 km altitude range. They observe 
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that In association with the warming phase of the stratwarm, there 
was found a reversal of zonal and meridional wind from westerly and 
southerly prevailing directions respectively. Gfiago^y and Mauon 
(1975) deduce meridional temperature gradients from their wind data 
and believe that the pattern of temperature changes In the vertical 
as suggested by Labitzke. (1972a, b; see next section) In relation to 
the stratwarms Is consistent with Saskatoon data. Other observational 
studies, for example by Ryazanova oX aJL. (1976), may also be cited 
for higher latitudes. 

In the present study, the radio meteor technique has been used 
to obtain wind data In the height range of 80-100 km over Atlanta 
(34 ®N, 84®W). The prevailing wind regime over this height range is 
analyzed for the period of August 1974-July 1975. Particular empha- 
sis Is placed on the circulation changes observed at meteor heights 
over this low mid-latitude station as related to the minor high lati- 
tude stratospheric warming that peaked during the first week in Jan- 
uary 1975. This warming event Is studied In the next section. A 
brief description of the meteor radar observation technique follows. 
Finally, the prevailing wind at the mesopause level is analyzed in 
subsequent sections. 

2. Stratospheric Warmings 
A. T/ie Mtno^ St^atoipkcJuc. Warning 1974-75 

Since its discovery by Scherhag in 1952, the anomalous warming 
of the winter stratosphere has come to be recognized as the single 
most energetic disturbance of the entire stratosphere. Individual 
warming events have been termed major or minor depending on whether 


large scale circulation reversals were observed In the middle strato- 
sphere (at or below 10-mb. ''30 km» level) In association with the 
temperature Increases observed (of at least 25° In a period of a week 
or less). In the higher latitudes {McJnttviUf 1978). 

The mid-winter warming of 1974-75 satisfied all but one criterion 
required to be classified as major ; that Is, the mean westerly cir- 
culation poleward of 60°N did not reverse to easterly at 10-mb level. 
Figures la and lb show 2-mb (''42 km) charts for December 11 , 1974 and 
January 1, 1975 respectively. On December 11, 1974 (Figure la), rep- 
resenting the pre-warming upper stratosphere, a cyclonic vortex is 
situated close to the pole and the cold center (dashed lines represent 
constant temperature contours) lies about 90° westward. The chart 
for January 1, 1975 (Figure lb) shows the state of the circulation 
close to the peaking of the warming event. Here the low pressure 
center Is seen displaced to a position east of Greenland at about 70° 
latitude and the Aleutian high centered over the Bearing Straits gov- 
erns the polar circulation. As shown in Figure 2, a plot of the mean 
zonal winds, easterlies are to be found north of 55°N at the 2-mb 
level on January 2, 1975. By January 15, weak westerlies appear again 
and within a week the stratospheric circulation returns to its near 
normal state. 

Figure 3, which is a plot of zonally averaged temperatures, shows 
that the warming reached its peak during the first week of January 1975 
in the upper stratosphere. On January 2, 1975, an increase in temper- 
ature of 40°C is observed at the 2-mb level, north of 70°N. Tempera- 
tures return to normal in the upper stratosphere only after the second 
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week of January 1975. Figure 3 also Indicates a low latitude warming 
In association with the cooling at high latitude (on January 8 and 15 
contours). VfUtz and Soules (1972) have shown that for the major 
warming of December 1969-January 1970 In high latitudes based upon 
satellite radiance data (weighting function In the middle stratosphere), 
a simultaneous decrease Is observed In the low latitudes of both the 
hemispheres and more Importantly, a phase lag Is observed In the radi- 
ance maxima from low to high altitude. 

Warming pulses of about 15-day periods have been found In the 
winter stratosphere of both hemispheres, and are readily detected 
from the satellite borne radiance measurements. In Figure 4 the mean 
zonal radiance (solid line) at 82.5°N Is plotted for the winter of 
1974-75. The peak marked II, during the first week of January coin- 
cides with the "minor" warming Jescribed above. In all, six peaks 
can be discerned over the period of three months. 

B. TmpoAatuAz Chang&i up to SO bn 

During a stratwarm the height of the stratopause, defined as the 
temperature maximum, can vary from 20 to 60 km altitude in high lati- 
tudes. The mesopause height however, seems to vary over a much nar- 
rower range, possibly from 70 to 90 km altitude. 

Lab^zke (1977a) has combined several soundings from different 
winters to give profiles for high and low latitudes for a "composite 
warming" (for early winter (I), climax (II) and late stage (III)l as 
shown in Figure 5. Here, opposite changes in temperature between strat- 
osphere and mesosphere (right part. Figure 5), and between high and 
low latitudes are apparent. The region of sign change lies between 
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50-55 km and Labitzke call it "the pivotal level for the interaction 
between stratosphere and mesosphere in winter". 

3. Wind Observations Using Meteor Radar 
A. The QQjoKQia. Tech MdAjo \slind VacAJLity 

The Georgia Tech system (at Atlanta; 34°N, 84®W) has been in rou- 
tine operation since July 1974. It has been designed as a continuous 
wave, all sky system with a capability of continuous operation 24 hours 
a day, seven days a week with an adequate usable echo rate (RopeA, 1975a, 
b). 

A double sideband suppressed carrier continuous wave transmitter 
operating on 32.5MHz+360Hz with an RMS output of 2 Kw has been installed 
on the Georgia Tech campus, and a receiving site established at Technology 
Park/Atlanta, 27 Km northeast of the campus. The Georgia Tech ranging sys- 
tem uses phase comparison of the sidebands of a double sideband suppressed 
carrier signal to determine the echo range. For direction finding, rel- 
ative phases between receiver outputs are determined by integration 
over 3/2 cycle of Doppler waveform, ensuring that each echo accepted 
meets the requirements placed on echo arrival angle. Data tapes at the 
receiving site contain information only on the Doppler frequency and 
relative phases, along with a continuously updated record of year, day 
of year, hour, minute and second. A series of FORTRAN computer pro- 
grams handles the data reduction from determination of echo arrival 
angle, range, trail drift velocity and time of occurrence of echo to 
the plotting of height/time wind profiles and altitude dependent wind 
spectra deduced from the echo data. These programs have been documented 
by RopeA (1975a) . 


B, VzcompoAition oi the. Meteo/i Wind 

The Georgia Tech facility uses a technique developed by G^ov&i 
(1957), that finds a "best fit" model to represent data by generalizing 
the least square solution. In this method, the zonal, meridional and 
vertical components of the model wind at the reflection point on the 
trail are assumed to be certain specified functions of height and time 
containing arbitrary parameters. These parameters are then chosen in 
such a way that the error between the model wind and the observed 
trail drift velocity is minimized. 

Typically, a polynomial variation in height, with a periodic var- 
iation in time, is allowed. Third degree polynomials in height for 
zonal and meridional components, and a first (or zero) degree polynom- 
ial for the vertical component are found satisfactory for data that 
has at least 120 echoes per grouping. The fundamental period is taken 
as 24 hours with 12 and 8 hour components also computed. This choice 
is justified by the ample data on the presence of large tidal compo- 
nents (24- and 12-hour) at meteor heights. 

The prevailing wind then, is the meteor wind averaged over a 
period of 5 to 10 days from which tidal components (24, 12 and 8 hour 
periods) have been removed. The echo rate dictates the period of aver- 
aging beyond five days. Figure 6 is the height-time plot of the zonal 
and meridional components of the prevailing wind for the period of 
August 1974 through August, 1975. There are no data available over the 
period of March 28-May 14, 1975. In this figure, isotachs are plotted 
at TO m/sec intervals. The continuous lines represent westerly and 
southerly winds (taken positive) on the zonal and meridional plots 
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respectively. The computed error (one standard deviation; not shown) 
over the entire period under consideration Is season dependent. It Is 
large during winter due to large variability of the wind at meteor 
heights during winter. It averages to about 20 per cent for a reason- 
able echo rate (I.e. about 100 usable echoes per day). The echoes 
are sampled over the height range of 76-106 km. though the winds below 
80 km and above 100 km are not considered due to large errors Involved 
(usable echoes outside 80-100 km height range are very small In num- 
ber). The height resolution Is 2 km In these experiments. 

4. Results and Discussion 
A. VkzmoXXXmq Wind, Auguit J974 through JuJLy 1975 

Zonal component . The zonal circulation Is characterized by mostly 
easterly wind during fall tfirough mid-winter. Westerlies are present 
during late winter and throughout the summer months (Figure 6). 

The slope of the zero Isotach during August 1974 Indicates a 
downward progression of the westerly regime. The westerlies reach the 
80 km level during the first week of September and appear to weaken 
while descending. A maximum westerly wind of more than 20 m/sec Is 
reached above 95 km around the third week In August. From mid-Septem- 
ber, easterlies re-establish beginning from lower levels, so that by 
the last week of October, the easterly wind takes over the entire 80- 
100 km height range. The easterlies gradually Intensify until a max- 
imum Is reached by the fourth week of November 1974 near 95 km; a peak 
value of greater than ,- 30 m/sec is recorded. 

Toward the end of November, easterlies start to weaken until a 
reversal begins In the third week of December 1974. The slope of the 
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zero Isotach, again* suggest a downward progression from 100 to 80 km. 

The easterly wind re-establishes, starting at lower heights by 
the beginning of January 1975, and intensifies to a peak value of more 
than -20 m/sec above 90 km during mid-January. 

Subsequent reversal to a westerly regime during the last week in 
January 1975 appears to have started from higher levels. Throughout 
February and March 1975, the zonal circulation is such that the west- 
erly wind strengthens with increasing height, while below 85 km, a 
weak and variable wind regime persists. 

There is no data for the period of April through mid-May 1975. 
From mid-May, where the 10 m/sec isotach is seen almost vertical, the 
wind is westerly through the summer and it strengthens with increasing 
height. 

The Meridional Component . Compared to the E-W component, the 
north-south wind is weak and variable in strength. 

The zero isotach located near 85 km divides the northerly wind, 
above, from the southerly wind, below, during August through mid- 
October 1974. During November, the zero isotach shows a progressive 
establishment of southerly wind over the entire 80-100 km height range 
starting at 80-km. The southerly wind persists through the last week 
in February 1975, with a few reversals to northerly wind. 

During March 1975, the wind is northerly and around the last 
week of March, there is a strong reversal (u > 10 m/sec) to southerly 
wind above 90 km. 

During the middle of May 1975, a strong southerly current is 
observed below 90 km (u > 20 m/sec) with northerly wind above. During 
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June, the meridional wind Is northerly, followed by southerly wind In 
July, From the third week of July, a reversal to wind from the north 
begins at higher altitudes. 

The Mean Meridional Temps:rature Gradients . Vertical shear In 
the geostrophic wind Implies the presence of a horizontal temperature 
gradient. The thermal wind, which Is the vector difference between 
the geostrophic wind at two pressure levels, cannot strictly be applied 
In relation to the meteor wind between the 80-100 km altitude range. 
However, qualitative Information about the mean meridional temperature 
gradients can be inferred from the height-time plot in Figure 6. 

The thermal wind blows parallel to the isotherms, with warm air 
to the right facing downstream In the Northern Hemisphere. Then, for 
example, when westerly wind gains strength or easterly wind weakens 
with Increasing altitude, positive thermal wind (zonal component) Is 
indicated; this implies that the mean meridional temperature gradient 
is negative and thus, colder latitudes are indicated nortliward. 

If at the same time, the mean meridional wind component of the prevail- 
ing wind is southerly, we can infer warm advection. 

Studying Figure 6 in this context, it can be seen that from mid- 
August 1974 through July 1975, the meridional temperature gradient is 
negative, implying cooler latitudes north of Atlanta, with two signif- 
icant reversals. During the period of the second half of November 
1974 and later on around mid-January 1975, warmer northern latitudes 
are indicated. 

From the last week in November 1974 through the last week in 
February 1975, the meridional wind is southerly at 80-100 km. Thus, 


during this period warm advection is indicated except around mid- 
January where there is cold advection at meteor heights. 

B, T/ie CVtcuZation Changes in the, UeJtzon. Zone Reiated to the 
SUatoAphe/Uc Ci^cutation SyAtm 

In this section the meridional and zonal components of the pre- 
vailing wind at meteor heights are analyzed in relation to the pressure 
and temperature regimes in the upper stratosphere. The period under 
consideration is from August 1974 through July 1975 (Figure 6). Major 
emphasis is on the minor stratwarm in the mid-winter of 1974-75. The 
effect of quasi -periodic warming pulses felt in the winter stratosphere 
is also studied. 

In Figure 7 are plotted temperature traces on Atlanta longitude 
(84 ®W) at 35®N and 50®N, at 0.4- and 5-mb levels, the temperature values 
being read from the synoptic charts produced by the Upper Air Branch of 
the National Meteorological Center 1978). The values are approx- 

imate and limited significance can be attached to the 35°N to 50°N tem- 
perature gradients read from this figure. It is also necessary to 
keep in mind that the thermal systems in the atmosphere generally change 
in intensity and position with changes in altitude. 

The Effect of a Strong Harming Pulse Felt During November of 1974 . 

The upper air synoptic charts 1978) indicate that during the 

second half of November 1974, a warming pulse was strongly felt in the 
upper stratosphere. Polar latitudes at 0.4-mb (~55 km) are warmer by 
about 35 °C on November 20, 1974 than on November 6, 1974. This can 
be seen from Figure 8, where temperature traces at 70°N latitude are 
plotted for November 6 and 20, and December 11, 1974 from 0.4-mb maps. 


A warm cell has moved over the Aslan land mss to a position 90°E 
within the Arctic Circle during November 1974. The 0.4-mb chart for 
November 20, 1974 (Figure 9} also shows two cold cells, one at 40°N 
off the China coast and another over Central Europe. There Is a low 
pressure center over Scandinavia and a weak low Is centered over the 
Northern North Pacific. At 5- and 2-mb levels (not shown), a weak 
high Is found centered over North Alaska ahd a weak warm cell can be 
detected westward at the 2-mb level (Sta^j}, 1978). 

Figure 7 shows a warming of about 10®C during the tflrd week 
of November 1974 at the 0.4-mb level for Atlanta's longitude. Also 
note the strong negative meridional temperature gradient (35®N to 
50°N) at the 5-mb level (~36 km). 

At meteor heights during the second half of November 1974 (Figure 
6), the easterlies strengthen above 90 km. The meridional temperature 
gradient Is positive, Implying warmer northern latitudes. The meri- 
dional wind is from the North above 92 km. It then appears that we 
have a thermal structure Indicated that will severely weaken the circu- 
lation associated with the stratospheric-lower mesospheric polar cyclone 
(I.e., westerly wind) that extends at least up to 100 km In stable win- 
ter conditions. 

It Is then probable that at middle and high latitudes, the warm- 
ing pulse that was felt strongly near stratopause heights extended up- 
ward, possibly beyond the mesopause (80-85 km). This conclusion fits 
well with Lab-Uzke's (1972a, b) principle of compensating changes of 
temperature in the vertical, when we note that below 40 km (Figure 7) 
there Is a tendency towards cooling in higher latitudes during the 
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second half of November 1974. 


The Effect of the Mid-Winter Warming of 1974-75 . The reversal 
to westerly wind toward the end of December 1974 and return of easter- 
lies in January 1975 leading to an intensification during the third 

I 

week, are the principal zonal circulation features that appear to be 
directly linked to the various phases of the mid-winter warming (1974- 

75) of the stratosphere. 

I 

We have reviewed this warming in Section 2. As mentioned before, 
this was a major warming event in the upper stratosphere and as Figure 

10 indicates, on January 1, 1975, at 5- and 2-mb levels (36 to 42 km) 

I 

the Siberian Arctic has warmed by about 40®C. Its maximum descent was 
probably around January 5, 1975 as the radiance peak in Figure 4 would 
indicate. 

The temperature traces for Atlanta longitude (Figure 7) toward 
the end of December 1974 show cooling at 35°N in the upper stratosphere. 

This compensating cooling in association with high latitude warming is 

I 

in line with the latitudinal compensation reported by VfUtz and Soules 
(1972). 

What are the implications of this amplification stage of the 

I 

stratwarm in relation to the circulation at meteor heights over Atlanta? 
From Figure 10 we notice that ihe warm stratopauso has descended to 
about 35 km height in polar latitudes by the first week of January 

» 

1975. This situation, in principle, corresponds to Stage II of Labitzke.'s 
(1972b) profile for "composite warming" (Figure 5), This profile leads 

to a very cold mesopause near 80 km height for higher latitudes; and 

I 

to a warmer mesopause in tropical latitudes (Figure 5; profiles on the 
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right). 


At meteor heights we find westerlies, strengthening with height, 
during this period (a reversal from easterly wind). The meridional 
wind is southerly and more intense during the second half of December 
1974 (Figure 6). The thermal wind is positive during December 1974, 
implying negative meridional temperature gradient over the entire 
height range, i.e., northern latitudes are cooler. 

Thus, it can be safely concluded that the circulation changes 
at meteor heights over Atlanta during the period of late December- 
early January of 1974-75 correspond to the changes implied for these 
heights by the warming State II described before. The meridional tem- 
perature gradient at meteor heights is conducive to westerly circula- 
tion. 

It is interesting to observe that the descent of westerlies from 
above 100 km during late December 1974 (Figure 6) appears comparable 
to the descent of easterlies (or warming; Figures 2 and 3) in the 
upper stratosphere, with the westerlies descending (implying gradual 
cooling) about a week (or more) earlier. 

Now taking up the late phase of the stratwarm, the zonally aver- 
aged temperature values (Figure 3) indicate that the warming of higher 
latitudes is considerably reduced by January 15, 1975 and the low lati 
tude warming is evident in the upper stratosphere. On Atlanta longi- 
tude, there is warming at 35°N and 50°N all through the 5- to 0.4-mb 
(36 to 55 km) height range. Also, by January 15, 1974, westerlies 
appear at all latitudes with a strong westerly jet centered near 40- 
45°N at stratopause heights (Figure 2). Figure 10 indicates warming 
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with Increasing height at 70®N, with a warm cell at 90®E. 

This situation, in principle, corresponds to Stage III (Figure 
5) of a major warming as projected by Lab-Uzke (1972a, b). As remarked 
earlier, the 1974-75 mid-winter warming was a majoi’ stratwarm in 
all but one regard— that the circulation reversal to easterly wind 
did not take place in higher latitudes at 10-mb or below. We then 
expect that at inesopause altitude, there will be a positive merid- 
ional temperature gradient, with a warm mesopause at higher latitudes- 
favorable to anticyclonic circulation at those heights. 

The observed changes in the zonal wind at meteor heights over 
Atlanta around mid-January 1975 (Figure 6) appear to support the pre- 
ceding remarks. The easterlies reappear in Janaury 1975, ascending 
from below (i.e., 80 km), and intensify during the third week above 
90 km (>20 m/sec). The thermal wind is negative, implying warmer 
northern latitudes. The meridional wind component has weakened in 
strength (southerly) through January 1975. 

The establishment of westerly wind towards the end of January 
1975 (Figure 6) from higher heights (i.e., from above 100 km) is the 
return of the normal circulation at meteor heights after the dissi- 
pation of the stratwarm. The cyclonic circulation in the upper strat- 
osphere is well established (notice the strong negative temperature 
gradient on Atlanta longitude. Figure 7). The cold mesopause has made 
a return in higher latitudes as implied by the strong positive zonal 
component of the thermal wind at 80-100 km. 

The Summer of 1975 . The stratospheric circulation during summer 
is governed by the relatively stable circumpolar anticyclone with east 


erly wind over much of the Northern Hemisphere. At meteor heights 
over Atlanta however, the zonal winds follow the thermal balance. 
Westerlies are present at 80-100 km (June through early August) with 
weaker winds at lower heights. Positive thermal wind implies a colder 
higher latitude mesopause. The meridional wind is weak and averages 
over the summer as weakly northerly as suggested by various models 
as well as other observations (for example, Uovy, 1964). 

C. GmoJwJL Comme.nt& on the. PAevcUtlng Wind StmxctuJie. 

It was observed earlier (section 4. A) that the prevailing wind 
(particularly the zonal component) generally exhibits a phase lag 
with increasing altitude. It is also observed (Figure 6) that the 
wind components are generally stronger at higher altitudes in the 
meteor zone. These observations suggest an upward propagation of wave 
energy from the atmosphere below into the region of decreasing density, 
at least up to 100 km. 

Next, taking a summary view of the implied meridional gradients 
(section 4. A), it appears that a cold mesopause is indicated through- 
out the year at northern latitudes. The only indications to the con- 
trary, during winter, have been related to the thermal perturbation 
of the high latitude stratosphere. In recent years Quifioz (1969), 
Labitzke (1972,b) and others have questioned the need of a heating 
requirement for the winter mesopause at high latitudes. Their specu- 
lation is that the winter mesopause may in fact be "cold" in early 
winter and again in the late winter, and the warm mesopause observed 
by rocket data (Th&on et at., 1972) is probably associated with the 
sudden warmings of the high latitude winter stratosphere. Present 




I 


i 


work Is in line with this interpretation. 

Finally, in this analysis the meteor wind has been interpreted 
in terms of the large scale circulation features of the stratospheric 
circulation system. In spite of the very local accumulation of the 
data, it is interesting to see how well the variation can be interpreted 
as part of the large scale circulation, and at such a relatively low 
latitude station at that. 

5. Conclusions 

In the present work, synoptic scale variations in the prevailing 
wind over Atlanta (34®N, 84“W) at 80-100 km are interpreted in the 
context of the stratospheric circulation system. The period covered 
is from August 1974 through July 1975. It has been shown that quasi- 
periodic warming pulses felt in the winter stratosphere can influence 
the wind regime in the meteor zone. The minor mid -winter stratwarm 
felt in the high latitude upper stratosphere during 1974-75 is studied 
and various phases of this warming have been shown to have affected 
the prevailing wind in the meteor zone over Atlanta in a consistent 
manner, by primarily making use of the latitudinal and vertical temper- 
ature compensation criteria in the stratosphere and above. What makes 
this analysis particularly interesting is the fact that the perturbed 
high latitude stratosphere interacts strongly with the mesopause level 
winds over a relatively low latitude station such as that at Atlanta 
(34®N). 

The large southerly wind during winter and weak northerly wind 
during summer at 80-100 km altitude is consistent with other observa- 
tions and mesospheric models (for example., by Leovt/, 1964). It also 
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appears that the Implied meridional temperature gradients In the one 
year data point to a year round cold mesopause at high latitudes. 

The only exceptions during winter are shown to be related to the per- 
turbed state of the stratospheric cyclone due to the warming episodes. 
Needless to add, data from more northward stations Is required to con- 
firm this Interpretation. 

Some wave energy Input from below Is Indicated by the prelimi- 
nary examination of the height-time variation of the prevailing wind 
structure. This Inference along with a possible explanation of the 
presence of easterlies during the fall of 1974 (where westerlies are 
expected to be present) are the topics of a future publication. Finally, 
It Is obvious from the present work that the radio meteor technique 
offers an Invaluable source of continuous data at lower themospheric 
heights for synoptic and long term studies of the general circulation 
system of the stratosphere-mesosphere. 
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Figure la: 2-mb chart for December 11, 1974. Isobars are solid 

lines at intervals of 320 meters, and isotherms are 
dashed lines at intervals of 5°C (after Staffs 1977) 
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Figure 2; Mean zonal wind associated with the mid-winter warming 
of 1974-75 {after Klinkey^, 1976). 
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Figure 3; Mean zonal temperature associated with the mid -winter 
warming of 1974-75 {af.ter Klinker, 1976). 






Figure 4: Rddy geopotential flux (GF) through the 100-mb “^surface 

integrated from 20“N to the North Pole (dashed line) and 
zonal mean radiance at 82.5®N based on VTPR (Vertical 
Temperature Profile Radiometer) Channel 2 measurements 
(solid line); 1974-75 (after Quiroz et at., 1975). 
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Figure 5: (a) Vertical temperature profiles at high latitudes for a 

composite mid-winter warming and also the changes between 
I and II and between II and III, (b) Vertical temperature 
profiles for Thumba, India. Dates as shown (after Labitzke, 
(1972b). 



70 


IAUG74 I SEPT I OCT I NOV I DEC ]JAN75 IFEB IMAR lAPB IMAY Uti. NE LIULY 


jliSAL PA-CfS IS 

pOOB. QUALITY 


Figure 6: Height -time variation of the prevailing wind components 

August 1974-August 1975. Wind speeds in m/sec. 
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Figure 9; 0.4-mb map for November 20, 1974. Explanation is in 

Figure la. 






temperature (®C) 



Figure 10: Temperature traces at 70“N latitude in the upper 
stratosphere on January 1 and 15, 1975. 




CHAPTER VIII 


MAP - The Middle Atmosphere Program 


\ 


i 




This report would not be complete without reference to the Middle 
Atmosphere Program of the Special Committee on Solar Terrestrial Physics 
(SCOSTEP) of the International Council of Scientific Unions. Details of 
this program, which Is designed to be to the middle atmosphere what GARP 
Is to the troposphere and the IMS Is to the magnetosphere, may be found 
In the MAP Planning Document (available from Dr. S. A. Bowhill, Department 
of Electrical Engineering, University of Illinois 61801). A recent (1980) 
update on the status of MAP In the U.S. Is available as "The Middle Atmo- 
sphere Program - Prospects for U.S. Participation" prepared by the Committee 
on Solar Terrestrial Research of the Geophysics Research Board, and published 
by the National Research Council, Washington, D.C. 20418. 

While the major thrust of MAP has been delayed until 1982-85, pre-MAP 
programs are underway. Pre-MAP 1, which has as Its focus a better under- 
standing of the effect on the upper atmosphere of winter polar stratospheric 
warmings, was held during the northern hemisphere winter of 1978-79. A con- 
certed effort to better understand the phenomenon of anomalous radio wave 
absorption In the D region is planned for 1980-81, with the West Germans 
spearheading an International rocket program to make In situ measurements 
In the mesosphere during the winter of 1980-81. Concurrent measurements of 
upper atmosphere dynamics using rocket released vapor trails and meteor 
wind radars have and will continue to provide vital knowledge of the dynamics 
of the upper mesosphere and lower thermosphere, a knowledge crucial to the 
understanding of these "anomalous" events. 
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Publications produced under Grants NsG-304>63 and NGL-1 1-002-004 


Evening Twilight Winds from 68 to 140 Kilometers for May 21, 1963, J. 
Res . . 1^, 6062-6063 (CiG. Justus, H.E. Edwards and D.C. Kurts). 


(Letter) 


Photographic Instrumentation for Triangulation Studies of Luminous Clouds In 
the Upper Atmosphere, Appl . Optics . 3^, 399-403, 1964 (M.M. Cooksey, Z. 
Frentress, H.D. Edwards) . 


AitSTKArr; The position ami velocity of clouils or oilier objects tiiovinji through 
the upper atmosphere are often ileterniineii from a sliniy of simultaneous photo- 
yrophs taken from two or more ohservinn stations anainst a ilarkeneil sky with a 
baik^roniul of siars. ll'ith eorreetion for atmospheric refraetion snpplieil by 
tahli the star harke,ronnil is ntilizeil to increase the aeenraey of camera orien- 
tation by correetinj!, errors inherent in the eamera system proiliieeil by: ' I ) him 
shrinlnijic, iiyht refraction in the yjass fuincial t>riii plate, and (,,<) tilt of thre 
e,tmera with resin et to the local horinonliil. Teehiibines and enipirical formnlae 
a e lie, eloped for use in analytic data procession with a diyitiil compnirr to the 
ordei of a thon.saiidth of a centimeter. The final procedure can produce anynlar 
position determinations of inilliradian, .1 method for accurately deter min- 
inn camera focal leitj>(l/s is also pre.sented. 


A Method Employing Star Backgrounds for Improving the Accuracy of the Loca- 
tions of Clouds or Objects In Space, Photoqrammetric Engr. , 2, 594-607, 
1964 (C.G. Justus, H.D. Edwards and R.N. Fuller). 


Instrumentation lias been developed for obtaining triangiilalion pliotograplis of luminous elouds in the 
upper ntmospliere. These eluiids are ereated li,v the release of ehemieals siieli as sodium and eesiiim 
from rockets. 'I’lie insirumenlation has liuen used liii appro.\imatel,v eighty rocket firings. Ueseriptions 
and operating charaeterislii's are given of the ei|iiipment. l tliservitlions have I'een made front the .Air 
Korce ninges at Kgliii .Air l''orce M.ase, l''lorid!t, and Ifolloman .Air l''orce Hii.se, N<‘\v Mexico, ami from the 
N.ASA lest facility iti Wallops Island, Virginiii. .A nioililied M'rsion of the iierial ciimcrii with 
7 iit. (17. S cm)/.’J..'i optii's was iisnl iiiid mounted on liO-in. ( l.'i'J-i'in') .searchlight carriages mid .Mk 'il 
gun dirt'ctors. (itinera control system, tlata chamber, tidiicial light system, shutter and (liter system, 
field installation, and alignnaait proeeilures are desi'ribed. 
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Ionospheric Winds: Motion Into Night and Sporadic E Correlations, Space Res, 
n/, 171-181, 1964 (N.W. Rosenberg, H.D. Edwards and J.W. Wright). 


Abstract; Tint puraistitncw uf night wiiula ami ahttura for more than H hoiira uftur siiiiRnt 
at ultitmk’R between lUO ami 150 kni wua ovatouteil (tiiring the COHVAK Byiio|tti<! 
ionuHpluii'ie wimla |ii'ogriim in December lOtbi. The ima,b(Ht iiseci wna to relimno 
ebeiniliiininiMm trails IVnio roekota between 100 ami 150 kni with the times of 
release being 1720, iMOl, 2115 ami 2245 OM'I' on K Deetanlter 1002. The release 
times wei')^ res|ieeliv’e|y HI minutes before Himset at the 100 km level amt 
2H minutes, 4 hours 12 minutes ami 6 hours 12 minutes pitst sunset. Alt roekets 
were launebed from tbe Air l'‘oree tiulf 'IV-st Kangit, Kglin .Air l'’oree Muse, Kloriila 
at a longilmie of HO.U" \\' ami a latitude of 20.0'' N. 

It is believitd that Ibis is the lirst time su<4i a juor|)hology of winds and sbeurs 
tbrougb till' /‘.'-region bus been made available. (Jomparisons nr" made between 
the wimis and shears observeil under night eomlitloiis with tbo.se for tbe preueisling 
twilight period. 

'I'beoretieal diseussk.ns by Whitehead r/ al. have suggested a wind shear origin 
for sporadie A,'. High east-west hori/.ontal wind shears in the ionosphere and the 
natural mairnetie tiekl will net on ioni/.ed speeies to stratify these into relatively 
high verlieal gradients. 

Tbe pre.sent paper reports on simultaneous measurements of ionospberie wind 
shears (measured by ebemilumineseenl. trails) and of sporadie I'l (by sweep- 
frei|ueney ionosomles). Deeiirrenee at the same altitudes of jiiasimum wind shears 
and strong sporadie A.' Ie\els was found. In one ease, for esampli> (SO.ti" W, 
20.ti’' N, l7‘wU t'ST, Deeeinber ,‘l, lOliw) sporadie A. was found at tMl {. I Uni lo 
I.U Me and at Itlh | I km at tl.2 Me. At this time, wind shear maxima of 
'15 m/s km bearing H5" K iit U7 ;l I km and of tio jn/s km bearing 27.5" K at 
109 .1; I lull were found. 

These data U>nd to verify the correlation between wind shear and sporadie K. 


A Triangulatlon Technique for Linear Objects In Space, Photogrammetrlc Engr. , 
31 , 1020-1029, 1965 (H.P. Haney, W.M. Schofield and W.H. Wooten. 


.\iisi i< act: Corresponding points on films, token from two difierent sites, of u 
linear object in space cannot be paircit doe lo unknown fores/iorloiint;. The. 
asnal Iriannidiilion techniques most be replaced by an iterative scheme to solve 
the problem. The particular methods developed in this paper are most useful 
when graphical checks of analysis 'ire desired. l''iirtliermore, the methods are. 
tailored for easy computer proy,ramini<^ and the closure technique of the ih ration 
allows easy and quick converyynce in pairing points. These methods have been 
used on chemical cloud releases in space to determine winds and wind shears. 
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Space and Time Correlations of Ionospheric Winds, Radio Sci 1, 149-155, 1966 
(C.G. Justus and N.U. Rosenberg). 


lonii8|ilu'rit' wiml |iuiii>riii» liHwi'vii 90 uimI 150 km wi'ri* Hltiiliml liy tiiNimiiitii o|' 
vapiM' truils Irom l'•ll‘kl’ll>, I'iltctMi Irails, I'litM'ly in kidk'i* or liini-. aiv ri’iioilml. l'ro> 

lili'.o .liotv »i>||-oijiaiii/.<‘<| .iiiiiMiiilal iiallmiK of iH‘ai'i‘om‘iaMl vvuvi-Ii'IikiIi. if llo’ vvrtii'al t'oonliiiali' 
Ih mfa’<(iM‘il in iol<‘i;rai<‘il mmIc lii’ittlih rallivr llian kiloioi'li'iMi Ikiiiiiiiaiil rvavo|i’m!;il> )<' - ‘f >^i'ali’ 
lo'iplil>. similar in NS ami l‘!\V i'om|HimmiM, ami pliasi' |iro)ir<'Ks*'s <|owiiivaril aliool I wavrlimatli (■•'r 
nir!lil, lli-ri/oiitally s|iai i‘il trails slio» vnlical |M'olili' ilis|ilai'mm'OI i’orri's|iimiliop to horiroiilal si ales 
ol UlOO 'I(KH) km ill NS wiml. loiiticr ill K\V (irolili-s. VisroiiH I’lii'i'tty ili»ni|ialioii irom oIisitvi'iI ivavos 
is a sittiiiltratil Imal soiirro at l.'iO km allilmli-, 


Energy Balance of Turbulence in the Upper Atmosphere, J, Geophys. Res. , Jl* 
3767-3773, 1966 (C.G. Justus). 


Turbnlont winds obtained from chemical release studies and turbulent tlilYtision observa- 
tions of globule expansion can bo used to olitain esUinntcs of llio terms iii llie I'ni'iay balmn’e 
cf|iialion e. = »» + ta, wtiere «. is the rate per unit imfss at wliicli enoriJiy is Mipiiln d by wind 
shears, and anil *j lire the rates per unit mass at wiiieli eiiurgy is rlissipated by buoyaiiey 
and viscous elTccls. Both «, nnd f, arc founil to vary slowly wilh nlliludi!, haviiig values of 
i\boul O.IS mid 0.1 wall k«"', respeelively, in llio idliludu rantn! 00 to 110 km. However, r.i in- 
eroa.ses rai'idly wilh nllilmle. This raind increase in ej is responsible for (ho nlobulu eiilolT 
observed at ail averaai! allilmle of 100 km. (Slolmln growth .mml.vsis indieales limt exUeiiH ly 
rapid anomalous dilTiisinn lakes place during the period upproximately 100 to 150 seconds after 
relea.se. As globule diameters increase boyoiul the largest scale of the anomalous growth region 
they exiiaml at a rate coniparablu to that of molecninr diffusion for a pcrioil of approximately 
SO seconds, after which a transition to d* I* turbulent diffusion occurs 


Dissipation of Wind Energy in the Height Range 80 to 140 Km, J. Geophys. 
Res., 71, 4427-4428, 1966 (R.G. Roper). 


(Letter) 
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Atmospheric Turbulence in the Meteor Region, J. Geophys. Res.. 71. 5758-5792. 
1966 (R.G. Roper). 


In 1001 at Adelaide tlie atmospheric tiirbuleni'e near 00 Inn was nicasiircd munth I'V month 
usiiii; a spacc<l-s!atiun radio meteor tcihiiiiitio, 'J'hc tiioriietcrii$lics of tliu hirac-M-iilo liirlndcnt 
rriuti(m.'r are found to be siiiiilm' to ihosm observed in the northern huinj.s|ilmio. (JlmraeteriMtii; 
rcloi'itk's range from 33 to <10 niMe and the large eddies are distinctly aniMitiopici they have 
a horizontal scide of from 00 to OSU km and u ycrtii'al scale of 7 km, in keeiimg with a gravity 
wave simrcij of tmlmli nt. eiiclgy. Slii’ais nii>ai>iM't d over M'|iaralioiiK of ((/» (o d.fi km are fouinl 
to be cliaraclerislh’ of an opiiarcntly inertial spoelruni, and llic tmlmhait, dissipation rate 
calenlaled from these shears shows a marked seasonal variation, with a toaviinuni of tUU 
ergs/g sec in autumn and spring and a minimum of 100 e.rgs/g sec in siimmi r and winter. 
There i, a strong corn lalion lietwoeii the seasonal variation of the lurli(d< at dissipation rate 
and the energy of iJio 31-hour component of the mean wind. The lu'ight sliear is of particidar 
interest, as, in keeping witli l•ocket vapor trail wind analysi.s, it follows aii ni)iiro.\imale l.d 
power law. 


The Eddy Diffusivities, Energy Balance Parameters, and Heati>q Rate of Upper 
Atmospheric Turbulence, J. Geophys. Res, . 7^, 1035-1039, 1967 (C.G. Justus) 


Wind proiile and turbulent wind dura obtained from photograpiric tracking of rocket re- 
leased ehcmieal clouds have been used to eompute the iliermal and momeutum eddy dif- 
fiisivities A'l and Km in the UO-ltO-km region. Ksponentiul functions having a small ioi'rea.ic 
wilit altitude provide a rcasoualric representation of the mean Ki and Km values. The ratio 
Km'Ki, the turbulent J'randll iminbei', is fuuial to be approximately constant with a value 
of about 3, Thu mean Ki at 100 km is uhout ISO iiiVsec, in leasomible agieciiniii with pre- 
viou.s e.slimales. The ex|iimeuUal fimction.s foi K, and K,„ are also used to olitaiii revinud 
e.slimales of the energy lialance parameters and the tmiioleiicu la atmg rate in tliis altitude 
range. 


The Spectrum and Scale of Upper Atmospheric Turbulence, J. Geophys, Res., 
72, 1933-1940, 1967 (C.G. Justus). 


Turbulent winds determined by plmtographie tracking of chemical release clouds are used 
to detui'inine tlie turbulence structure function in tlie UO-llO-km region. Tlie observations 
indicate that the turbulence structure function is approximately isotropic with a large scale 
of about 5 km, 'J’lie turbulence strueturi' fmii'tiun is found to vary, :is r"' in tin 1- to .'i-km scale 
range. However, there is insullieieiil. observational data um| theorelieal luicKgroimd to delor- 
miuu if I Ills is a I rue inertial scale range or is a Tsemlo-incri ial’ region in wliiih wind sliear 
produelion ami biioyiomy lo.ss are approximalely bnlam'cd, leading to no m l ht.ss from llm 
S|)e<‘trum. .Structure fum'lions of llii) total wind proiile are also discussed and related to a 
pu.ssil)le form for the giavity wave s'lieetrum. 
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The Velocity Probability Density of Upper Atmospheric Turbulence, J. Geophvs. 
ReSj_, 72, 2460-2462, 1967 (C.G. Justus and W.B, Moseley), 


(Letter) 


Turbulence in the Lower Thermosphere, with M. G. Elford, Space Res. VJI. 
42-54, 1967 (R.G. Roper). 


Abstract. Results of studies of wind shears in the lower thermosphere (SO- 120 knt) 
by means of observ.ntions of ohemiluminoscent trails and meloor trails are 
summarized. The motion is highly anisotropic - in (he heiuht ranKc 80-100 km 
typical scales are 2<)0 km horizontal and .‘>-10 km vei’tical. The horizontal mo- 
tion, s cxliildt the characteristics of an inertial subrnnt;c of eddies. 

The rale of dissipation of turhulent cncri?y increases cx|H)nentially over llic 
heiuht r.anKC 80-lo.i km. with an average value of 3 x 10"® VV/kg at !)3 km. The 
efiuivalont heating rates of tlie atn»osphere over this same hciglu range arg 1-.»®K 
per day. Turbulence decreases r.api(lly above 10.> km. 

Below luu km the meteor observations indicate that there is a strong seasonal 
variation in, the rate of dissipation of turbulent energy; at 03 km a minimum value 
of 1.5 X 10"^ W/kg occurs in summer and winter, and a maximum value of 
3.5 X 10'^ W/kg during spring and autumn. 


Reply (to Discussion of a Paper by C.G. Justus "Energy Balance of Turbulence 
in the Upper Atmosphere"), J. Goophys. Re s.. 73, 45'j-4bH, 1968 (C.G. Justus). 


(Letter) 


Atmospheric Tides in the Height Region 90 to 120 Kilometers, J. Geophys. Res ., 
73, 467-478, 1968 (C.G. Justus and A. Woodrum). 


Winds determiued from chciiiical release trails in the UO- to 120-kin region (primarily from 
the early und middle fall season) are analyzed by a method timt is an exteii.^idii of tlic 
method useil by Ilincs fur tlio delcrminaliun (if the diurnal tide. Tlii.s niellmd is ile,<iaimil to 
reveal the prevailing, diiirnnl, and scniidiiirnnl componenl.s. A (linnial tid(> is compnu d iliuL 
has a vertical wavelength (HI ± 1 km) in agrceinem with tim meaMiremenls of llims and 
with tidal tlieoiy preiliei ions fur (he (1,3) mode, 'riii' observed lieialii. viiiialimi of ilie pre- 
vailing wind can mo.st ea.sily bo e.vpiained in tornis of a wave whose vertical wau'lenglli ia 
27 .-i: 5 km. 'J'lie t.inall observed magnitude of the seniiiliiirnal tide inilieales that there is 
eoiisiderable seasonal VLirialioli of this cunipotieiii, since a largo seniiiliiirnal tide has pre( ioii.sly 
been observed diiiiiig perimls near the .snninier mid winter solstices. Observed pha.-o dili'ei- 
(•ncos between the nortliward ami e:isl\vari| eompoiients of llie prevailing and tidal wimis are 
diseiissod. Some sugge.stions for iinpruvements in the tidal exlrirelion analysis method are 
a bo given. 


Some Observations of Turbulence in the 90 to 110 Km Region of the Upper 
Atmosphere, AMS Meteoro. Mono. . 31,, 122-128, 1968 (C.G. Justus and R.G. 
Roper). 


AHHTUAt r 

Some rwfiu mt'.iHiircmeiiis of Imbulcmu in llte 80-110 km rf(jion oln.iiiiwl by both the r.nlio-metuoi ami 
I'hemical release mc(ho(U are »iimmarizeil, The velority |)rob.ibility density, the sirneture fimetioii, eddy 
dilftisivities, heatiOK rate, ami the altitude ami seasonal v.irlation of the eiierny balance par.imeters of the 
iiirbulenrc are disrusHcd. Some as|)ects of the aiiparent disaKreeinent lieineen meteor trail ami chemical re- 
lease measarenienis are prescuied, and a pro|>osed ipialitaiive model is olTereil for the behavior of iiirbaleace 
in this reKion of the atmospliere and for the iiuerprelalion of the experimental results. This moilel ascribes 
die chemic.d release wimi lluctnalions ami the snl.ill-sc.ile meteor wind dillerences to liiibnlence, but at- 
tributes the larm'-scale meteor wind lliiciualions to gravity wave .iction. |)issip.itibn of tidal and itr,ivity 
wave enerny is (uoposed ;is a sonire of the turbulence eiieiyy, timl some .ispects of the w.ivi'-to-hu biilem e 
iransfer .ire diM'n«sed. 


Measurement of the Magnitude of the Irregular Winds in the Altitude Region 
Near 100 Kilometers, J. Geophys. Res, . 73, 7535-7537, 1968 (C.G. Justus 
and A. Woodrum). 


(Letter) 


Measurements of Tidal Period Winds in the 95 to 135 Km Region, J, Geophy s. 
Z 5(J , 6 , ) . 4099-4104, 1969 (C.G. Justus, A. Woodrum and R.G, Roper). 


A mellioil developed by O. V. (iroves for tlie niial.t .sis of inelror wiml.s Im.s bi'en ii.'cd 
oil ehetniral release wind dal.-i Irttin roi'kels. Veloeily aiiiiditiide.s. [phases-. ,-,i)d ain)>liiiaie erroi-.s 

for winds willi peiioils of ,S. 12. and 21 I were eoin|mh'd in Ha- !l,a- to la.a-km reuioli. An 

avenme wind, a.ssuined |o reiu'esi’iiL a ronstanl |ll■l■^lolln^i: wiml was al.-o epiin|pn|p-p| in lids 
alliliido I'anne. Veilii-al wavelenylh.s and pha.sn dilTp'i’ein-e.s of the norihwai-p| and i-aslwaid 

l■om|muenl.s of Ihn iiri.-vailiiiK and tidal winds ale e.-deidali'd and l•oln|lin■ell wilh l■e.snll.s. 

l)l-.si|ialiM,' foree.s iire tilso shown to be aelinj; on the projiauiptini’ dinmal and sp niiplini'ind 
tides in the heitild. range from Ha to 1115 kin. 


Dissipation and Diffusion by Turbulence and Irregular Winds Near 100 Km, 
J. Atmos. Sci. . 26(5) , 1137-1141, 1969 (C.G. Justus). 


AUSl'RACr 

A review is j?ivcn of some recent measuremems of the (allowing energy lialimee imramcterB in the Wi-1 10 
Itni region: the viscous dissipation of liinetic energy by slieurs in the mean wimlsj the viscous 

rlissipaii.m of kinetic eiuTg>' by shears in tlie lurlHilent wimlsj Itie transfer of kim iic energy from lliu 
mean motion to the tmlmleni motioiij llie transfer of lurbiilenl kinetic energy to poienlial energy l>y 
Inmyaney action; ami the rale of ilissipalioii of wave energy of tides ami irregular winds tnterpreleil as 
gravity waves. Some measurements of the growili rales of globular Hlnictiire on r liemical release tiomis 
and the growth rale of inlerglolailar distances are presented. The.su data Indicate that tin: dilfnsion nieeli.i- 
nisin is a niislure of the variaiioo e.\pected sor the dilTiision of a point from its initial position with the 
,.(/* variation espeeled for the variation of the .separation between pairs of points. 


A Theory for the Energy Spectrum of Shear- Dependent Turbulence, J. Atmos. 
Sci., 26(6) , 1238-1244, 1969 (C.G. Justus). 


AHSTK.Srr 

,\ llieorv tor llic th.'ce diioensiiiiial energy .spectriim of oearli isolrojiii slu ,ti-di pcMl. m iiirluili la i is 
presented. Tliis theory is based on tlie author's gi'iieral moilel for shear tli'tiemlent Itirlmlem e .ind a moditied 
I’no spectral transfer iheorj wliieh acconnis for the elTeels of viseoiis loss, shear prodm lion, ,ind im rli.il iind 
velocity graiiii’iit transfer in nonslalionary luriniicnce. Tlie lestillanl spectral iransfer l•'|(|,|(ion is s.ilvcd 
in cluseil form for the cast^ .d cmniileie similarity. 'I'he small wavemiinl.ir st.ilionary u ergy -pce'.iiiin is 
simwn to vary as Tlie small wavemindier minstalionarv solnlion varies as wlu're t) ran liavi v.dms 
from 1-4 Oiu'-ilimensioiial s|ieelra are computed mmnrically ami compared witli tiie Uirlmleine spectra 
nliserved for wimi tunnel Inrlmlenec wilii a wide range of Reynolds mimlwrs. The theory successfully 
lirediris several features of the observed spectra at both iiigli .nml low wavemimbcrs as well as in the inertial 
region, 


Distribution and Structure of Irregular Winds Near 100 Km, J. Geophys. Res. . 
75, 2171-2178, 1970 (C.G. Justus). 

I’ho veioeify dialribntiun and strneturc function of irregular winds near tlic IflO-kin level 
aru ijrnsnntnd. 'J’ho iiregnlar winds wori! inefiiiired liy ililTereneing r'lminieal ridtm.sn wind 
proftli.'.s dolermined at limes ililTerent, tiy tin inleginl iniilliplo of 24 lioin.s tind by eoinimting 
Uio deviiilions from tlie 3-da,v mean of winds at llm .same lime of d.'iy. There is n signilieant 
vaiiaiion of llm rm.s irregular winds la lweeii the !M)- and ll.'i-km limil.s over whieli lh<y were 
nie.'i.'-iiml. The iiiagmlmle.s .and v.iii.ali(in.s willi altiliido me fomid to he rea.->mialily .-imil.ar 
for dal.i ohiaiimd I'roii, Kglin AI M, I'lorid.i; Viima, Arizona: .iiid llarli.ulos, ^Ve^t Indh'.s. The 
dislrihiitions u( irregid.ar wiiids seei.i lo <Je|iart slightly fioni a t.atis'.iim tlistrilmlion. In 
|iarlieidar, it low value of tlie llatness ftiHor lias been measiiied eonsialenllv. 'I'lie idliinde timl 
lime >lfne.tmv fuiielioiiH of the in'egnlaf winds indieate a elmraeleristie. verlie.il wavelength of 
alnnit Id km and lime st ales of ii)i lo 20 limns. 'I’he obserted slnielnre fmielions over small 
si aifs (2-.'l kio for (lie vmlieal sirmdnie fmielion, aboiil 2 hours for the lime slmetn.e fme 
lion) follow a power law variation with an expoiieiil- value near Jl/2. Tiiis may re.siiJt frtmi ;m 
energy .speelrnm of irregular wiml.s (Imt varies as the —f>/2 power of n;ive mimla r and 






Response of Winds in the 90 to 140 Km Altitude Region to Variations in Solar 
Activity, J. Geophys. Res. , 75^, 5565-5570, 1970 (C.G. Justus and J.E. Hicks). 


Winds df’lprminrd from S2 |)hotO)jrn|iliiciil|.v •ninkcd IriiilH of roi'kol-liorm' ilHnoK'dl ri'loiiscs 
ill till' !)0- III HO-kni idliliulc rnnifi', Iniinrliod lii'Uvi’i'ii Oclolicr of l!Ki2 mid Jiiiii' of 
iiri' iisi'd ill lioth n ((fiTi liilioii miidysis mid ;i liiionr I'l'Kriwioii miidysis iviili hi vcrid wdiir mid 

Iji oiiiii|»iM'lir iiidiri's, Ui'MiIIn indii’iiti' llini the iiKiliidc iiilmiil loiihiili ird iii.iy 

Ih' lirolo'ii n|i into two dislini'lly ililTi'iont ii'Kioiih of siiliir iidlii'’i|i‘i' ; llio icuioiih iilmvi’ mid 
lioliHV 110 I.III. Sinnilii’mil iiM'iviiKCH in tlio ui'oln.' ii'lii' iiid< N A'/» ovir ii iri-lmiir intiTvid 
l■^’llll•l'l'll on II linii' liiK of 21 lionr.s nfli'r iiii'ii'iiscs n 11)1110, siihcnV rircidiOion liolmv 110 km 
wlKiohl I'illmr 11 dynmiio ly|H' of iiili'rmOion or mi iniiinsi' in llic moiiiiiyimlir fiidd 
nIi'i'iikiIi, inodiii'i'd hy I'lirimsi'iiliir rudiiilion wliii'li lii(j.s lii'liind tlio l'll'l•lronlll^!nl■lil• rndiiiiion 
mid iimini|imi,viiiK ii(iiios|dii'iii' I'iri'iilnlioii viirinlioiis, Ai'i'inviitly ronlinnoii.s iissoi'iiOion.s, 
ivilli jiidiciilod iii'i'iod.s of iilioiil oiii' sidnr loliitioii, liidwoi’ii Inilli llio sidiir rndio iuiisi< index 
/'' 111.1 mid till’ Ziirii'li .siiiiHiiol niinilx'r Kt mnl tin* norlliwiird wind I'oniiioiii'iil inimiiiludr iiimve 
no km indii'iiii' tliii) loiiB-lcrin viiriiilions in tlie winds in lids iiBion me l•^■llll(‘d lo similnr 
vmimiiiiis in soliir KUV riidiiiliiin. 


Thermospheric Observations Combining Chemical Seeding and Ground Based Tech- 
ri^quest, Planet. Space Sci . . 20 , 761-789, 1972 (D. Rees, K.H. Lloyd. C.H. 
LOW, B.J. McAvaney, aiid R.G. Roper). 


Alixlriicl- Two Skylark sounding rockuls carrying chemical seeding payloads wure launched 
fniiii Woomera, Soiilh Australia (3TS, I37"t;) in October 1969. In ciinjunction with these 
lii'iiigs, the University of Adelaide eondueted gnnmd-based e.\|ierinieiits on the upper atmos- 
phere using the radio meteor end spaced receiver drift methods. This paper pre.seiits the nieus- 
III ements of properties of the neutral atinosphcre above 90 km which were obtained from these 
e.speriineiils. 


Upper Atmospheric Planetary Wave and Gravity Wave Observations, J. Atmos. 
Sci ■ , 30 , 1267-1275, 1973 (C.G. Justus and A. Woodrum). 


AbSTK.VCT 

Previously collected data on atmospheric pressure, density, tcniiieralure mid \i inds lictween 2.’ and 200 km 
from smirics imludini,' .Meleoroloaiial Umket N'elwork data, UOItlN’ falliiiK sphere darn, prci adc iclcase 
and pilot lulic data, meteor «imls, i hemieid release sviiids, saiellile data, and olliers acre analt.'.ed In a 
daily dilTerem e mellmd, and resulls on ilie mat;nilude of atmo.splierie pertiirlialioos inleipreled as eraviiy 
uaves and planetary waves are presented. Traveling plmielary-wave iiinliilmiions in the 2,s S.s km raiiue 
were foiinil to have a siBiiiliiant Ih'IbIh and lalitiidinal varlalinn. It was fouml llial olisi-rved eraviu wave 
density perturliainms and winil are related lo one anotlier in llie manner pre<lii led by eravily-wave ilieon 
ll was delernimed llial, on the averaae, uravily-wave energy deposilion or lelleelion oi i nrs at all .illiliides 
estepi llie .35 -75 km region of the iiiesospln're. 


The Measurement of Meteor Winds Over Atlanta, Radio Sci.. 10. 363-369, 1975 
(R.6. Roper). 


An ull-iiky, contimiuuh wave radio meteor wind facility has been installed in Atlanta by the 
Georgia Institute of Technology under National Science i’oundatjon sponstHship. A double sideband 
suppressed carrier CW tiansmitkr. operating on ^2,5 Mllr. ' Will II/, with an rius output of 
2 kw, has been installed tm the (ieoigia Tech campus, and a receiving site establi.chcd at Technology 
Park/ Atlanta, 27 km northeast of the campus. Details of the etiuipment, bigether with height /time 
piofiles of mean wind cireulutiou and tides between KO tind KMI km, measured from -.ugusi 'i 
to September 4, l‘>74, are presented. 


A Global Reference Atmospheric Model for Surface to Orbital Altitudes, J. 
Appl . Met e orol . . 15 , 3-9, 1976 (C.G. Justus, A. Woodrum, R.G. Roper, and 
O.E. Smithyr 


AltSTRACT 

An empirical atmost>hcric moilel has been developed which generates values fur pressure, density, Icm- 
perultire and winds front surface levels to orbital altitudes. The output parameters consist of components for : 
I) latitude, louKitiide, iiml altitude dependent niontlily means', 2) i|uasihientiial oscillations; and .1) ranilom 
perturbations to partially simulate the vuiiabiiity dite to synoptic, diurnal, planetary wave and gravity wii.ve 
vuriutiiiiis. The niontlily lOeuii models consist of : (i) NASA's four dimensional worldwide model, tieveloped 
by ICnvironniental Kesearcli and Tetluinl>igy, for height, latitude, and longitude ilepcndeiit monthly means 
from the surface to 25 km; and (ii) a newly developeit latitude-longitude dependent model which is an ex- 
tension of the Groves latitude depetideiit model for the region lielween 2.5 and IM) km. The Jucchia IV70 
model is used above CO km and is faired with the inodifieil Groves values between 'Ml and IIS km. Quasi- 
biennial and runilom variation perturbations are computed from purunieters determined from various empiri- 
cal studies, ami are added to the monthly mean values. This model has been developed as a computer program 
which can be used to generate altittttle profiles of atmospheric varialdes for any month ut any desired loca- 
tion, or to evaluate atmospheric parameters along any simulated trajectory through the atmosphere. Various 
applications of the model are discussed, and results are presented which show that good simulation of the 
thermodynamic and circulation characteristics of the atmosphere can he achieved with the model. 


Turbulence in the Lower Thermosphere, Chapter 7 of "The Upper Atmosphere and 
Magnetosphere", F.S. Johnson, ed., one of the Studies in Geophysics pub. 
by the National Academy of Sciences, Washington, D.C., 1977 (R.G. Roper). 


(See Chapter 2) 
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Winds from the Atlanta (34“N» 84“W) Radio Meteor Facility, J. Atmos. Terres t. 
Phys. , 40, 891-894, 1978 (R.G. Roper). 


T 


Absiricl— A brief description of the Georgia Tech radio meteor wind facility is followed by a tabular 
prcsuniaiion and discussion of winds measured over Atlanta (34°N, K4°W) for the first three intervals 
of the URSI/IAOA Cooperative Tidal Observations Program (CTOP), The pervailing zonal wind 
nieasured during August P>74, being easterly, is signilicantly different from that measured during 
Ocioher l*t7S and January l')76, and is not typical of winds measured in August l*>7S and August 
Id7f>, when westerlies predominated. The complicated tidal picture is detailed, but is nut easily 
summarized. 


A Comparison Between Radio Meteor and Airglow Winds, J. Geoniag. Geoelectr. . 
31, 419-426, 1979 (G. Hernandez and R.G. Roper). 


A comparison between the winds near 97 kilometers altitude has been made 
from observations of the I7924K (5577A)()( line emission at l•lll/. Peak Observatory 
(39,fsN, |95.s\V)aiul with a meteor radar facility at Atlanta (.14N, S4Wj, fiom An* 
gust 1074 to November l‘>75. Since the optical emission measurements are made 
only at night, the nighttime meteor ra.lar measurements hate been used, weighted 
by an airglow emission rate prolile. The results show general tigreenKiit in liotli 
the zonal and meridional wind vectors, btit with the variations in the amplittide of 
the meridional winds at the northernmost station (I'rit/ Peak Observatory) larger 
than those at Atlanta, a result of the smoothing inherent nt producing the meteor 
winds. 
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